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I.  INTRODUCTION 


A.  BACKGROUND 

Friis  and  Feldman  [1]  proposed  the  idea  of  an  electronic  beamformer  in  1937, 
when  they  published  “Multiple  unit  steerable  antenna  for  short  wave  reception,”  which, 
according  to  [2],  “supported  a  very  sophisticated  test  of  angle  and  time  diversity 
combining.”  Although  their  studies  led  to  the  development  of  modern  technologies  like 
multiple  beam  arrays,  digital  beamforming,  and  space-time-adaptive  processing  (STAP), 
they  were  limited  to  existing  technology  in  the  1930s  and  1940s  and  had  to  use 
mechanically  rotated  phase  shifters  to  develop  multiple  unit  steerable  antenna.  In  the 
1950s,  the  ferrite  phase  shifter  was  introduced,  which  led  to  improvements  that  allowed 
better  beam  steering  control  [2], 

In  the  following  decades,  radar  systems  took  advantage  of  electronically  steered 
phased  arrays.  Aperture  antennas  (e.g.,  reflectors)  are  efficient  and  simple  to  construct 
but  are  limited  by  mechanical  beam  control.  For  electronic  warfare  (EW)  applications, 
arrays  play  a  major  role  in  direction  finding  and  tracking  systems.  One  important 
example  of  a  phased  array  antenna  application  is  the  U.S.  Navy  AN/SPY-1,  shown  in 
Figure  1.  The  electronically  controlled  beam  steering  is  faster,  more  precise  and  versatile, 
and  allows  the  radar  to  search  for  and  track  multiple  targets  simultaneously. 

Phased  array  antennas  are  not  exclusively  employed  in  military  applications.  For 
example,  distributed  antenna  arrays  have  been  used  for  decades  in  long  baseline 
interferometry  in  astronomy  [3].  In  meteorology  applications,  weather  radars  take 
advantage  of  electronic  beamforming  to  improve  weather  radar  performance,  such  as  in 
the  MWR-05XP  mobile  phased  array  weather  radar  [4],  Recently,  automotive  collision 
avoidance  radars  using  phased  arrays  have  been  installed  in  automobiles  for  “automatic 
cruise  control,  parking  aid,  blind  spot  detection,  and  side  collision  warning 
capabilities”  [5]. 
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Figure  1 .  DDG-74  (USS  McFaul)  with  AN/SPY-1  Radar  Highlighted 

Source  [6]:  “Ballistic  Missile  Defense:  The  Aegis  SPY-1  Radar,”  Mostly  Missile 
Defense,  2012.  [Online].  Available:  http://mostlymissiledefense.com/2012/08/03/ 
ballistic-missile-defense-the-aegis-spy- 1  -radar-august-3-20 1 2/.  [Accessed:  23-Nov- 
2015], 


More  recently,  the  Defense  Advanced  Research  Project  Agency  (DARPA) 
sponsored  a  program  to  develop  new  phased  array  hardware  technology  which  will  form 
the  basis  for  new  DOD  projects  employing  arrays.  The  Arrays  at  Commercial  Timescales 
(ACT)  will  research  and  expand  existing  technologies  to  develop: 

1.  A  common  building  block  for  RF  arrays, 

2.  A  reconfigurable  electromagnetic  interface  (the  antenna  interface  from  the 
electronics  to  the  RF  signal),  and 

3.  Over-the-air  coherent  array  aggregation  [7], 

If  successful,  the  DARPA  ACT  program  will  provide  a  starting  point  for  all  array 
development  programs.  The  main  factor  leading  to  expensive  array  projects  is  the 
engineering  costs  for  developing  a  new  product,  which  are  usually  specific  for  an 
application  and  platform  [7].  ACT  acknowledges  that  a  “technical  solution  is  needed  to 
bring  military  array  programs  to  more  manageable  cost  levels  and  timescales,”  [7]  and 
this  program  will  facilitate  the  access  of  this  technology  by  industry  to  manufacture 
civilian  systems. 
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B.  DISTRIBUTED  DIGITAL  PHASED  ARRAY 


To  achieve  a  narrow  pencil  beamwidth  and  high  gain  and,  consequently,  long- 
range  detection,  systems  like  the  AN/SPY- 1  are  physically  large  and  structurally  difficult 
to  handle.  The  system,  once  implemented,  is  not  flexible  (i.e.,  it  is  not  suitable  for 
significant  changes). 

To  perform  beam  steering  and  side  lobe  suppression,  the  array  beamformer  needs 
to  accurately  control  the  element  excitation  phase  and  magnitude  (i.e.,  the  element 
weights).  Note  that  the  following  discussion  is  with  regard  to  a  transmitting  array,  but  the 
same  concepts  apply  to  receiving  arrays. 

Traditionally,  analog  beamforming  requires  a  large  and  complex  network  of 
transmission  lines,  phase  shifters,  and  attenuators  to  control  the  radiated  waveform  at 
each  individual  element.  A  photograph  of  an  S-band  ground-based  radar  phased  array  is 
shown  in  Figure  2,  which  clearly  illustrates  the  weight  and  size  disadvantages  of  an 
analog  beamforming  network. 


Figure  2.  A  Phased  Array  Antenna 

Source  [8]:  “Phased  Array  Antenna,”  Hughes  Aircraft  Co.,  Culver  City,  CA  (unpublished 
photo). 
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Digital  arrays  use  a  digital  beamformer  (DBF)  controller,  which  is  a  computer, 
instead  of  the  analog  beamforming  network.  The  analog  portion  of  the  array  resides  in  the 
transmit-receive  (T/R)  module,  which  is  located  at  the  element  level.  The  differences 
among  these  architectures  are  illustrated  in  Figure  3(a)  and  3(b). 


(c)  Wirelessly  networked  digital  phased  array 


Figure  3.  Phased  Array  Architectures 

Source  [3]:  D.  Jenn,  Y.  Loke,  T.  C.  H.  Matthew,  Y.  E.  Choon,  O.  C.  Siang,  and  Y.  S. 
Yam,  “Distributed  Phased  Arrays  and  Wireless  Beamforming  Networks,”  Int.  J.  Distrib. 
Sens.  Networks,  vol.  5,  no.  4,  2009.  p.  284. 


In  a  distributed  architecture,  digital  arrays  can  have  their  elements  spread  along 
large  surfaces  like  the  entire  ship  hull  or  a  hillside  [3],  overcoming  some  of  the  problems 
that  arise  with  large  arrays. 
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In  similar  fashion,  the  wireless  distributed  phased  array  proposed  in  [3]  has  the 
beamformer  network  replaced  by  wireless  network.  “This  architecture  brings  a  set  of 
advantages  including  flexibility,  scalability,  adaptability  to  the  operational  environment, 
and  survivability”  [3]. 

The  DBF  connects  to  T/R  modules  through  a  wireless  local  area  network 
(WLAN)  using  commercial  off-the-shelf  (COTS)  equipment.  The  wirelessly  controlled 
phased  array  is  illustrated  in  Figure  3  (c).  The  wirelessly  networked,  digital  distributed 
phased  array  is  aligned  with  the  concepts  in  the  third  technological  area  of  ACT,  i.e.,  the 
development  an  over-the-air  coherent  array  aggregation. 

The  great  advantage  for  ship-based  digital  distributed  arrays  is  the  possibility  of 
connecting  subarrays  that  are  spatially  separated  through  a  wireless  or  wired  network; 
therefore,  the  possibility  of  synthesizing  large  arrays  on  a  single  platform  exists  using  this 
distributed  digital  subarray  (DDSA)  concept  [9]. 

C.  PROBLEM  DESCRIPTION 

Phased  arrays  are  designed  to  have  high  gain,  accurate  beam  pointing  and  low 
sidelobes;  however,  due  to  the  use  of  non-ideal  components,  magnitude  and  phase  errors 
occur,  altering  the  pattern  characteristics  if  these  errors  are  large  enough. 

There  are  many  sources  of  phase  errors,  and  all  components  of  the  system 
contribute,  direct  and  indirectly,  to  increase  these  errors.  Timing  and  phase  errors  related 
to  local  oscillators  (LOs)  are  the  main  source  of  phase  noise  and  frequency  drift,  and 
LOs  are  used  as  a  reference  for  clocking  in  digital  systems;  hence,  phase  noise  and 
frequency  drift  affect  all  digital  components,  such  as  direct  digital  synthesizers  (DDSs), 
analog-to-digital  (A/D)  converters  and  digital-to-analog  (D/A)  converters. 

There  are  other  sources  of  phase  errors  that  are  not  directly  related  to  LOs. 
Quadrature  mixers,  which  usually  are  integrated  into  modulator  and  demodulator  boards, 
can  distort  phase  information  if  an  exact  90°  phase  shift  between  in-phase  (/)  and 
quadrature  (0  LO  is  not  supplied.  For  example,  homodyne  modulators  use  independent 
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DDS  channels  to  generate  I  and  Q;  consequently,  a  difference  in  D/A  gain  can  impact  the 
phase  balance.  The  same  occurs  in  homodyne  receivers  [10]. 

Moreover,  phased  array  antennas  rely  on  the  coherent  processing  of  transmitted 
and  received  signals  among  their  elements  for  beamforming;  therefore,  phase  errors  are 
introduced  if  the  element  position  is  not  correctly  processed. 

The  objective  of  this  thesis  research  is  to  find  acceptable  limits  on  amplitude  and 
phase  errors  in  the  up  and  downconversion  processes  by  examining  their  effects  on  gain 
loss,  beam  pointing,  pattern  distortions,  and  degradation  in  sidelobe  levels  based  on 
parametric  simulations. 

In  most  instances  it  is  assumed  that  all  errors  are  uncorrelated.  It  is  important  to 
mention  that  there  are  several  types  of  highly  correlated  errors  which  can  increase 
sidelobe  peaks  significantly.  The  frequency-dependent  phase  errors  due  to  transmission 
lines,  quantization  errors,  and  imprecise  tapering  may  generate  frequency  and  spatially 
correlated  errors  [11]. 

D.  RELATED  WORK 

For  the  past  several  years,  many  students  have  contributed  to  distributed  digital 
phased  array  research  at  Naval  Postgraduate  School  (NPS).  One  of  the  most  critical 
challenges  to  overcome  is  the  problem  of  time  and  frequency  synchronization.  Loke  [12] 
and  Tsai  [13]  investigated  wireless  local  oscillator  distribution  and  element 
synchronization  using  a  “brute  force”  algorithm.  In  [13]  Tsai  achieved  phase 
synchronization  with  a  20-degree  tolerance  range. 

Additionally,  Loke  [12]  also  investigated  the  impact  of  deforming  surfaces,  i.e., 
ship  hull  deflection,  in  studying  element  mislocation  effects  on  synchronization.  For  large 
wavelengths  such  as  the  VHF/UHF  frequency  bands,  it  was  found  that  no  correction  was 
needed;  however,  at  higher  frequencies  phase  errors  due  to  position  errors  can  be 
significant. 

In  2008,  Djerf  and  Tomazakis  [14]  designed  and  tested  T/R  modules  for  a 
distributed  digital  array  radar  (DDAR)  and  also  developed  a  control  and  processing 


6 


program  in  Lab  VIEW  to  control  elements  using  a  field  programmable  gate  array  (FPGA). 
They  also  addressed  phase  synchronization  problems  by  investigating  different  T/R 
module  components,  such  as  circulators,  power  dividers  and  low  noise  amplifiers 
(LNAs).  In  their  experiments  synchronization  was  achieved  with  10  degrees  of  deviation. 

In  addition  to  phase  error  introduced  by  the  synchronization  algorithm,  addressed 
in  [12],  [13]  and  [14],  other  aspects  of  phase  synchronization  such  as  jitter,  frequency 
drift  and  frequency  offset  contribute  to  distort  beam  pattern  and  sidelobe  levels.  These 
phenomena  are  described  in  [15],  where  a  mathematical  model  for  instantaneous  phase  is 
described. 

Array  error  effects  due  to  random  magnitude  and  phase  errors  were  analyzed  in 
[11]  by  Mailloux.  Analytical  expressions  were  derived  to  evaluate  the  impact  on  residual 
sidelobe  errors,  peak  sidelobe  expectation,  gain  degradation  and  beam  pointing  error. 

E.  THESIS  OUTLINE 

This  thesis  is  arranged  into  six  chapters.  In  Chapter  I,  general  background  about 
phased  array  antennas  is  presented.  The  thesis  objective  and  problem  description  were 
also  provided,  along  with  a  summary  of  previous  work. 

In  Chapter  II,  the  T/R  module  architectures  and  the  issue  of  phase  errors  in  the 
frequency  multiplication  stage  caused  by  LOs  errors  are  presented.  Additionally,  in  this 
chapter  is  summarized  the  research  on  LO  distribution  and  synchronization  techniques 
published  in  [12],  [13],  [14]  and  [16], 

In  Chapter  III,  digital  array  theory  is  reviewed.  The  expressions  for  array  pattern, 
array  factor  and  gain  are  presented,  along  with  a  modified  array  factor  model  that 
includes  phase  errors. 

In  Chapter  IV,  a  mathematical  model  for  non-ideal  LO  instantaneous  time  and 
phase,  which  includes  the  impact  of  phase  noise  and  frequency  deviations  is  presented. 
The  phase  noise  phenomenon  is  addressed  in  this  chapter  by  discussing  the  Power-Law 
and  the  Leeson’s  models. 
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In  Chapter  V,  the  simulation  results  are  presented,  and  the  effects  of  phase  errors 
on  pattern  distortions,  degradation  in  sidelobe  levels,  gain  loss  and  beam  pointing  are 
illustrated.  Additionally,  a  typical  scenario  is  analyzed  to  exemplify  the  effects  of 
amplitude  and  phase  errors  in  the  up  and  downconversion  processes  using  the  simulation 
results  presented  in  this  chapter. 

In  Chapter  VI,  the  research  is  summarized  and  conclusions  are  presented,  along 
with  suggestions  for  possible  future  research. 
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II.  CRITICAL  TECHNOLOGY  AND  TECHNIQUES 


A.  T/R  MODULES 

1.  Transmitting  Architectures 

Transmitters  can  be  divided  into  homodyne  and  heterodyne  types.  Homodyne 
transmitters  have  only  one  upconversion  stage  and  for  this  reason  are  also  called  direct 
upconversion  transmitters.  Heterodyne  transmitters  generally  have  two  upconversion 
stages,  which  take  baseband  signals  to  an  intermediate  frequency  (IF)  and  then  use 
another  mixer  to  convert  IF  to  the  carrier  frequency  at  a  radio  frequency  (RF).  The 
architectures  are  illustrated  in  Figure  4  and  Figure  5,  respectively.  The  digital  portion  of 
the  receiver  is  highlighted  with  a  dotted  box.  The  array  beam  control  is  performed  with  a 
digital  signal  processor  (DSP).  A  power  amplifier  (PA)  is  used  to  provide  the  necessary 
gain  for  transmission,  whereas,  a  band  pass  filter  (BPF)  is  employed  to  suppress  spurious 
signals. 


Figure  4.  Homodyne  Transmitter  Architecture 

Adapted  from  [10]:  L.  B.  Oliveira,  J.  R.  Fernandes,  I.  M.  Filanovsky,  C.  J.  M.  Verhoeven, 
and  M.  M.  Silva,  Analysis  and  Design  of  Quadrature  Oscillators,  vol.  1.  Springer 
Science+Business  Media  B.V.,  2008.  p.  16. 
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Figure  5.  Heterodyne  Transmitter  Architecture 

Adapted  from  [10]:  L.  B.  Oliveira,  J.  R.  Fernandes,  I.  M.  Filanovsky,  C.  J.  M.  Verhoeven, 

and  M.  M.  Silva,  Analysis  and  Design  of  Quadrature  Oscillators,  vol.  1.  Springer 

Science+Business  Media  B.V.,  2008.  p.  16. 

Complex  analog  baseband  signals,  i.e.,  I  and  Q,  are  generated  in  homodyne 
transmitters  by  a  DDS  in  two  different  branches.  In  the  modulator,  I  and  Q  signals  are 
mixed  with  the  LO  independently  and  then  combined  in  the  RF  front  end.  This 
architecture  requires  two  D/A  converters  and  a  high  power  BPF. 

The  main  constraint  with  this  architecture  is  the  local  oscillator  “injection  pulling” 
and  “injection  locking”  effects  due  to  poor  isolation  between  the  mixer  and  power 
amplifier.  Both  effects  cause  a  frequency  disturbance  in  the  transmitted  signal  [10], 

The  heterodyne  transmitters  combine  I  and  Q  signals  at  an  intermediate  frequency 
and  then  upconvert  them  to  RF  using  a  single  mixer.  By  using  two  stages  of 
upconversion,  heterodyne  transmitters  have  to  deal  with  image  frequency  rejection  and 
also  demand  high  attenuation  band  pass  filters  (around  60  dB)  to  reduce  unwanted 
spurious  frequencies  [10]. 
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2. 


Receiver  Architectures 


Receiver  architectures  are  divided  into  three  groups:  homodyne,  heterodyne  and 

low  IF. 

Homodyne  and  heterodyne  receivers  have  one  and  two  downconversion  stages, 
respectively.  Homodyne  receivers  downconvert  RF  signals  using  two  mixers  to  generate 
baseband  I  and  Q,  which  are  sampled  and  digitalized,  as  shown  in  Figure  6.  Heterodyne 
receivers  first  downconvert  from  RF  to  IF  and  then  use  another  stage  to  obtain  baseband  I 
and  Q,  as  shown  in  Figure  7.  A  LNA  is  used  to  improve  the  received  signal  power. 

Homodyne  receivers  are  sensitive  to  flicker  noise,  channel  selection,  LO  leakage, 
quadrature  errors  (/  and  Q  amplitude  and  phase  errors),  DC  offsets,  and  intermodulation 
(IM)  effects  [10],  and  consequently  might  not  be  the  best  architecture  for  T/R  modules. 


Figure  6.  Direct  Downconversion  Receiver  Architecture 

Source  [10]:  L.  B.  Oliveira,  J.  R.  Fernandes,  I.  M.  Filanovsky,  C.  J.  M.  Verhoeven,  and 
M.  M.  Silva,  Analysis  and  Design  of  Quadrature  Oscillators,  vol.  1.  Springer 
Science+Business  Media  B.V.,  2008.  p.  10. 


The  drawback  for  heterodyne  receivers  is  the  necessity  of  two  downconversion 
stages,  which  complicates  chipset  integration  and  requires  an  RF  filtering  stage  to 
suppress  the  possible  image  frequency  before  mixing  [10];  furthermore,  the  additional 
downconversion  stage  increases  phase  noise  and  spurious  responses. 
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Figure  7.  Heterodyne  Receiver  Architecture 

Source  [10]:  L.  B.  Oliveira,  J.  R.  Fernandes,  I.  M.  Filanovsky,  C.  J.  M.  Verhoeven,  and 
M.  M.  Silva,  Analysis  and  Design  of  Quadrature  Oscillators,  vol.  1.  Springer 
Science+Business  Media  B.V.,  2008.  p.  9. 

Combining  previous  techniques,  we  obtain  a  third  approach  with  the  benefits  of 
both  heterodyne  and  homodyne  receivers.  A  low  IF  receiver  is  a  heterodyne  receiver  that 
suppresses  the  image  frequency  in  the  mixing  stage.  An  example  is  the  Hartley 
architecture,  illustrated  in  Figure  8,  in  which  the  RF  signal  is  first  do wncon verted  to 
complex  IF  (i.e.,  I  and  Q ).  The  image  frequency  is  suppressed  by  adding  I  shifted  by  90° 
with  Q.  The  output  of  a  Hartley  mixer  is  a  clean  IF  output,  which  is  treated  like  regular  IF 
in  heterodyne  receivers.  This  architecture  requires  a  very  accurate  quadrature  signal; 
therefore,  in  the  case  of  1  and  Q  unbalance,  the  image  frequency  rejection  is  incomplete. 


Figure  8.  Low  IF  Receiver  Architecture 

Source  [10]:  L.  B.  Oliveira,  J.  R.  Fernandes,  I.  M.  Filanovsky,  C.  J.  M.  Verhoeven,  and 
M.  M.  Silva,  Analysis  and  Design  of  Quadrature  Oscillators,  vol.  1.  Springer 
Science+Business  Media  B.V.,  2008.  p.  12. 
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3.  Transmitter  and  Receiver  Components 

a.  Local  Oscillators 

Oscillators  are  used  to  supply  an  accurate  reference  signal  to  quadrature  mixers 
for  converting  the  RF  signals  to  baseband  or  IF  signals,  and  vice  versa.  Their  main 
characteristics  are  a  high  Q  factor,  low  phase  noise,  and  tunability  [10].  Beamforming 
literature  usually  deals  with  an  ideal  LO;  however,  an  inaccurate  LO  reference  causes 
distortion  in  the  beamforming  process  that  affects  sidelobes,  gain  and  beam  pointing 
angle. 

Oscillators  are  divided  in  two  groups  [10]: 

1.  Strongly  non-linear  or  relaxation  oscillators  are  usually  designed  as 
resistor-capacitor  (RC)  active  circuits,  a  feature  that  allows  them  to  be 
integrated  into  a  chipset.  They  are  often  used  in  applications  that  are 
resilient  to  phase  noise. 

2.  Quasi-linear  oscillators  are  those,  which  have  a  resonator  element,  such  as 
dielectric  resonators,  crystals,  striplines,  and  inductor-capacitor  (LC) 
tanks.  Resonators  give  the  circuit  better  performance  when  compared  to 
relaxation  oscillators. 

One  common  example  of  quasi-linear  oscillator  is  the  voltage  controlled  crystal 
oscillators  (VCXO)  using  a  quartz  resonator.  A  quartz  resonator  is  a  capacitor  using  a 
piece  of  quartz  crystal  as  dielectric.  Since  quartz  is  a  piezoelectric  crystal,  when  excited  it 
converts  electric  potential  into  mechanical  oscillation  [17].  The  advantages  of  a  VCXO 
are  high  Q  and  low  phase  noise,  which  is  desired  for  distributed  arrays;  however,  the 
frequency  selection  range  is  very  limited. 

b.  Mixers  and  Digital  Multiplication 

Mixers  are  non-linear  circuits  used  to  achieve  a  multiplication  of  two  signals.  In 
receivers,  mixers  downconvert  RF  signals  by  combining  them  with  a  local  oscillator 
reference.  Most  of  the  applications  use  balanced  mixers,  which  have  inherent  isolation  of 
RF,  IF  and  LO  ports  and  the  ability  to  reject  spurious  frequencies.  On  the  other  hand, 
balanced  mixers  require  a  large  number  of  components,  i.e.,  diodes,  and  require  high  LO 
power,  limiting  their  use  to  applications  without  LO  power  constraints. 
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Furthermore,  there  are  other  interactions  between  a  LO  and  mixer  that  must  be 
mentioned.  The  LO  amplitude  modulated  noise  is  added  to  LO  signals  when  it  is 
amplified  and  frequency  multiplied.  Additionally,  temperature  variations  and  vibration 
cause  noise  insertion.  These  random  errors  are  represented  by  LO  phase  noise,  which  is 
proportional  to  the  inverse  square  of  resonator  Q  factor  [18];  therefore,  phase  noise  is 
reduced  using  high  Q  oscillators. 

Another  concern  is  LO  spurious  frequencies  (or  spurs)  generated  in  mixing.  Spurs 
are  the  result  of  non-linear  combining  in  the  mixing  process.  Unlike  noise  spurs,  they  are 
usually  phase  modulated  and  difficult  to  filter  when  they  lie  in  the  RF  or  IF  bands. 

Intermodulation  (IM)  products  are  the  consequence  of  the  interaction  between  two 
tones  in  a  non-linear  device.  In  the  mixer,  the  two-tone  interaction  (i.e.,  LO  and  RF 
signals)  is  desired  to  up  or  downconvert;  however,  mixing  also  generates  other  unwanted 
product  frequencies.  For  example  the  output  in  the  mixing  process  gives  the  following 
frequencies  [18]: 


fo[F  =  mfo  +  nfw 


(2.1) 


where 

fQ  =  an  intermediate  frequency, 

fQ  -  radio  frequency 

flQ  =  local  oscillator  frequency, 

/ n,n  =  0,±1,±2...  oo. 

Additionally,  the  LO  spur  frequencies  f^Q  in  turn  generate  products  in 


spur 


do  wncon  version : 


(2.2) 


and  likewise  in  upconversion: 


forr=nfott,+’fu)r.- 


(2.3) 
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The  frequencies  f0m  and  fn  must  be  sufficiently  far  from  the  passband  to  be 

111  spur  spur 

efficiently  filtered  out.  For  example,  as  illustrated  in  Figure  9,  the  desired  signal  is  the 
lower  sideband  m,n  =  -1  which  crosses  through  the  passband.  No  other  products  enter  the 
passband,  and  all  are  sufficiently  removed  from  the  passband  to  be  filtered  effectively. 


Figure  9.  Spurious-Response  for  a  30  GHz  Mixer.  LO  Frequency  is  26  GHz. 

The  IF  Band  is  from  3  to  4  GHz. 

Adapted  from  [18]:  S.  A.  Maas,  “System  Considerations,”  in  Microwave  Mixers,  Boston: 
Artech  House,  1993,  p.  182. 


Third-order  intermodulation  should  be  avoided  for  heterodyne  receivers  since  co¬ 
channel  intermodulation  can  cause  interference,  whereas  the  second  order 
intermodulation  is  worse  for  homodyne  receivers  because  co-channel  intermodulation 
might  cause  DC  value  insertion.  A  balanced  mixer  combined  with  a  high  Q  oscillator 
might  avoid  some  of  these  issues  [10]. 

In  [10]  the  authors  present  a  concept  of  a  balanced  mixer  integrated  with 
oscillator  that  provides  very  accurate  quadrature  conversion,  which  reduces  undesirable 
intermodulation  products  and  additionally  provides  a  precise  LO  reference  for  I  and  Q. 
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c.  Direct  Digital  Synthesizer 

The  DDS  is  a  component  that  uses  digital  signal  processing  to  generate  analog 
signals.  DDS  digital  characteristics  allow  for  a  production  of  a  variety  of  signals  with 
precision  and  stability  when  compared  with  analog  techniques.  The  DDS  architecture, 
exemplified  in  Figure  10,  comprises  of  frequency  and  phase  accumulators,  a  lookup  table 
and  D/A  converters  [19]. 


Figure  10.  General  DDS  Architecture 
Source  [19]:  M.  I.  Skolnik,  Radar  Handbook,  3rd  ed.,  vol.  1.  McGraw-Hill,  2008,  p.  6.48. 


B.  SYNCHRONIZATION:  LO  DISTRIBUTION  APPROACHES 

A  digital  array  requires  accurate  phase  synchronization  among  all  elements  in 
order  to  achieve  high  performance  beamforming.  An  inaccurate  LO  reference  corrupts 
the  waveform  by  introducing  magnitude  and  phase  errors,  which  causes  distortions  of  the 
array  pattern. 

In  order  to  maintain  the  array  synchronization,  three  different  LO  distribution 
approaches  can  be  taken: 

1.  Direct  distribution  of  a  single  LO 

2.  Use  of  independent  stable  LOs  (STALOs) 

3.  Phase  estimation  through  a  data  network  clock 
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Direct  distribution  of  a  LO  through  a  wired  or  wireless  channel  is  the  simplest 
architecture  to  keep  remote  nodes  synchronized.  In  a  direct  distribution  scenario,  the  LO 
signal  can  be  a  single  tone  radio  frequency  transmitted  from  a  controller  to  the  elements. 
In  an  ideal  hardwire  distribution  scenario,  all  cables  should  have  the  same  length  in  order 
to  guarantee  the  same  phase  reference  for  all  nodes;  however,  this  implementation  may 
be  impractical  for  large  distributed  arrays,  demanding  several  of  meters  lossy  coaxial 
cables  installed  along  the  platform. 

Loke  [12]  and  Tsai  [13]  tested  a  wireless  LO  distribution  technique  using  variable 
phase  shifters  on  T/R  modules  which  consisted  of  increased  phase  delay  on  shorter  path 
elements  to  match  the  delay  of  longer  ones.  This  LO  distribution  technique  required  high 
power  LO  signal  transmission  to  compensate  for  propagation  spreading  loss  and  deal 
with  different  path  distances;  for  instance,  Tsai  in  [13]  used  the  modulator  board  AD8346 
that  requires  LO  input  power  within  -6  and  -12  dBm. 

In  2006,  Loke  [12]  investigated  a  “brute  force”  synchronization  approach.  A 
reference  continuous  wave  (CW)  signal  was  sent  by  the  digital  beamformer  where  one 
element  was  selected  as  master  and  other  as  slaves.  The  received  CW  from  the  master  is 
transmitted  back  to  beamformer  and  is  used  as  a  reference  signal.  The  synchronization  is 
performed  one  slave  element  at  a  time.  The  slave  element  gradually  phase  shifts  its  LO 
signal  and  transmits  back  to  the  beamformer.  The  beamformer  verifies  element 
synchronization  by  subtracting  the  reference  (master)  and  element  returns  and  searching 
for  a  null.  The  master  and  slave  are  synchronized  when  the  null  is  achieved.  If  the 
algorithm  is  performed  often,  then  slow  channel  changes  will  not  affect  array 
synchronization.  During  experiments,  Loke  [12]  showed  that  the  whole  process  requires 
only  a  few  iterations  to  achieve  sufficient  phase  synchronization. 

The  plot  in  Figure  1 1  is  the  residual  phase  error  for  simulation  of  synchronization 
for  an  array  of  100  elements  that  are  randomly  located.  After  convergence  the  residual 
phase  error  is  within  ±  1 1.25  degrees  using  four  bits  of  phase  shift,  which  represents  22.5 
degrees  per  iteration.  An  advantage  of  this  approach,  which  uses  a  single  common  LO,  is 
that  frequency  drift  is  common  to  all  elements. 
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Figure  11.  Residual  Phase  Error  After  Synchronization  with  4-bit  Phase  Steps 

Source:  [12]  Y.  Loke,  “Sensor  synchronization,  geolocation  and  wireless 
communication  in  a  shipboard  opportunistic  array,”  M.S.  thesis,  Department  of  Electrical 
and  Computer  Engineering,  Naval  Postgraduate  School,  2006. 


In  2008,  Djerf  and  Tomazakis  [14]  looked  at  two  techniques  for  achieving  stable 
local  oscillators  in  the  T/R  module:  (1)  using  phased  locked  loops  (PLLs)  or  (2)  using  use 
demodulator  boards  to  recover  I  and  Q  and  estimate  the  phase  through  I  and  Q 
magnitude.  A  PLL  has  the  ability  to  control  a  voltage  controlled  oscillator  (VCO)  based 
on  a  reference  input  signal,  which  can  be  a  LO  signal.  The  PLL  can  complement  the 
direct  distributed  LO  scheme  because  it  does  not  require  constant  LO  transmission  and 
has  the  possibility  of  transmitting  a  lower  frequency.  A  PLL  in  locked  mode  can 
synchronize  an  element  using  hardwire  or  wireless  link  connected  to  a  digital 
beamformer  using  the  process  explained  previously.  After  synchronization  the  PLL  is  set 
to  hold-over  mode. 

Djerf  and  Tornazakis  [14]  investigated  the  use  of  demodulator  boards  to  recover  I 
and  Q,  where  an  independent  LO  is  used  for  each  element.  For  ideal  components  and  the 
LO  operating  at  nominal  frequency,  the  baseband  I  and  Q  results  in  two  DC  valued 
signals  that  give  the  phase 


0  =  tan 


(2.4) 
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This  technique  is  sensitive  to  I  and  Q  unbalance  and  AM  noise;  furthermore,  the  LO 
operating  at  exact  nominal  frequency  is  unlikely.  In  this  case  the  phase  error  propagates 
to  other  components  of  T/R  module,  disrupting  the  array  pattern. 

Recently,  many  time  and  phase  estimation  techniques  for  wireless  sensor 
networks  (WSNs)  have  been  studied  using  phase  estimating  algorithms  through  time 
stamps  on  data  packets.  One  consolidated  time-estimation  protocol  used  for 
communications  to  data  stream  synchronization  is  the  IEEE  1588  Precision  Time 
Protocol  (PTP);  however,  the  accuracy  of  the  recovered  clock  is  still  not  sufficient  for 
radars  that  employ  frequencies  higher  than  VHF.  For  example,  in  [16]  Texas  Instruments 
reported  a  test  regarding  clock  synchronization  using  PTP.  For  the  experiment  a  slave 
DP83640  device  generated  a  10-MHz  nominal  clock  with  226  ps  mean  and  2.7  ns 
standard  deviation  using  information  transmitted  through  an  Ethernet  network.  In  other 
words,  for  a  10-MHz  clock,  95.4%  of  the  time  the  phase  deviation  is  less  than  20°  and 
clock  frequency  error  is  around  2200  parts  per  million  (ppm). 

C.  SUMMARY 

In  this  chapter,  the  basic  T/R  module  architecture  was  presented  along  with  the 
technological  challenges  for  digital  phased  array  development.  In  the  next  chapter  the 
digital  array  theory  is  reviewed  and  a  modified  array  factor  model  including  phase  errors 
is  proposed. 
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III.  DIGITAL  ARRAY  THEORY 


The  concept  and  theory  behind  digital  arrays  is  described  in  this  chapter.  The 
basic  expression  for  array  factor  is  used  for  both  analog  and  digital  arrays.  For  the  latter, 
signals  are  synthesized  at  each  element,  emulating  the  effect  of  a  beamforming  network, 
and  for  this  reason,  the  sources  for  phase  errors  are  different  than  for  analog  arrays.  A 
modified  array  factor  taking  into  account  the  architecture  specific  magnitude  and  phase 
errors  is  derived,  and  expressions  for  gain  and  directivity  are  given. 

A.  ARRAY  EQUATIONS 
1.  Array  Pattern 

The  total  pattern  array  is  the  joint  contribution  of  all  elements  of  the  array.  Each 
element  is  weighted  to  point  the  mainlobes  of  the  transmission  and  reception  patterns  in 
the  same  direction.  To  compute  the  pattern  it  is  necessary  to  specify  the  element  locations 
and  orientations.  The  frequency  domain  (phasor)  expression  for  the  pattern  factor  for  an 
array  of  size  N,  for  which  elements  are  located  at  positions  (. xn,yn,zn ),  n  =  1,  2,  3,. . .,  N , 
as  derived  in  [3]  and  [13]  is 


N 


F(d,$)=^WneJk-rnEFn(0,<p ) 


(3.1) 


where 


W 

" n 


weight  coefficient. 


xVn  =  VSn+Wn 


magnitude  of  complex  weight  for  the  element. 


Vn 


phase  of  complex  weight  for  the  element. 


beam  scanning  phase. 
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uv 


U 


W 


A 


scan  x  direction  cosine, 
sin  0S  cos^y 

scan  y  direction  cosine, 
sin  0S  sin  <j)s 

scan  z  direction  cosine, 
cos05 

scanning  angle, 

position  vector  from  array  origin  to  «th  element, 

[v+yny+znz] 

observation  x  direction  cosine, 
sin^cos^ 

observation  y  direction  cosine, 
sinflsin^ 

observation  z  direction  cosine, 
cos6> 

wavelength. 


The  complex  weight  Wn  is  set  for  beam  scanning,  sidelobe  control  and  phase  corrections. 

2.  Element  Factor 

The  characteristics  of  an  array  beam  pattern  are  also  dependent  on  the  electric 
field  radiation  pattern  of  the  individual  elements,  called  the  element  factor  (EF).  The  EF 
is  the  electric  field  pattern  due  to  the  current  flowing  on  the  antenna.  We  assume  that  the 
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array  is  large  so  that  mutual  coupling  variations  can  be  ignored  and  all  element  patterns 
are  identical. 

The  element  factor  of  a  general  cosine  type  pattern  in  the  direction  r  is  expressed 

as  [13] 


EF  =  \ 


'V  I  *  -aI  C( 

p\n»r\  , 
0 


A  A  A 

n»r>  0 
otherwise 


(3.2) 


where  the  exponent  q  controls  the  half-power  beamwidth  (HPBW).  The  unit  vector  p  is 
added  to  indicate  polarization  direction.  For  example,  for  an  element  beam  pointing  along 
the  --axis  (h  =  z),  a  general  far-field  linearly  polarized  element  factor  can  be  expressed 
as 


EF  =  EFp  =  EF  (6  cos  Q  +  Qsiruj))- —  (33) 

where 

EF  =  element  factor,  and 

p  =  polarization  vector. 

For  a  distributed  array  the  element  pointing  directions  and  polarizations  can  vary. 
For  a  planar  array  with  uniform  element  spacing  and  alignment,  the  polarization  vectors 
and  normal  vectors  are  the  same  for  all  elements.  When  calculating  the  array  radiation 
properties  only  the  normalized  EF  is  needed.  A  convenient  expression  for  the  element 
factor  is  a  cosine  to  a  power  q  as  defined  in  (3.2). 

Another  concern  regarding  element  pattern  is  the  field-of-view  (FOV),  which  is 
generally  considered  to  be  the  element’s  HPBW.  For  uniform  arrays,  all  elements  have 
the  same  FOV.  In  contrast,  arrays  distributed  along  curved  surfaces  have  elements  that 
might  not  be  able  to  illuminate  the  target  due  to  obstruction  or  limitation  of  their  FOV 
[3];  hence,  these  obstructed  elements  do  not  contribute  to  the  array  pattern  and  EF  =  0. 

Given  identical  element  factors,  (3.1)  can  be  reduced  to 
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F(0,Q)  =  EF 


f  N  -  -  3 

I  K‘Jhr" 

V«=l  J 


EFAF 


(3.4) 


where  the  sum  in  parentheses  is  the  array  factor  ( AF) 

N 

AF  ='LWneJk'rn  ■  (3-5) 

n- 1 


For  the  remainder  of  this  study  we  focus  on  the  properties  of  the  array  factor  and  suppress 
EF.  With  regard  to  (3.4)  we  consider  F  to  be  a  complex  scalar  equal  to  AF. 

3.  Array  of  Subarrays 

The  use  of  subarrays  is  motivated  by  the  difficulty  in  allocating  sufficient  area  of 
a  platform  surface  for  a  large  array.  Breaking  the  array  into  smaller  subarrays  and  having 
them  installed  in  noncontiguous  (separated)  areas  minimizes  the  impact  of  the  array  on 
the  platform  structure.  On  the  other  hand,  the  received  and  transmitted  signals  from 
multiple  subarrays  must  be  coherently  processed. 

Consider  a  collection  of  M  arrays,  m =1,2,3. .  .M  where  the  subarray  centers  are  at 
[. xs{ni),ys(ni),zs(ni)\ ,  as  illustrated  in  Figure  12.  Let  be  the  array  factor  for  the 

mth  subarray.  In  this  case  the  total  array  factor  of  the  combined  subarrays  is  expressed  as 
[9] 

M 

mf)  =  Z  Ftn  (0,f*)exp  JyA:([x5(mX«  -«*)]  +  0»Xv-  Y*)]  +  [zs{m){w  -  ws  )])}(3.6) 

m= 1 

where  Fm(9,<f> )  is  the  array  factor  for  the  subarray  m  as  given  in  (3.5). 
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Figure  12.  Subarray  m  and  its  Local  Coordinate  System  Relative  to  the  Global 

Origin  0. 

Source:  B.  Feng,  D.  C.  Jenn,  and  S.  Member,  “Two-Way  Pattern  Grating  Lobe  Control 
for  Distributed  Digital  Subarray  Antennas,”  IEEE  Trans.  On  Antennas  Propag .,  vol.  63, 
no.  10,  pp.  4375-4383,  2015. 

4.  Including  Phase  Errors  in  the  Array  Factor 

The  precision  of  pattern  control  depends  on  two  main  factors:  the  application  of 
the  correct  value  for  complex  scanning  weights  Wn  for  each  element  and  the  knowledge 
of  the  exact  location  of  the  elements;  therefore,  the  use  of  non-ideal  devices  that  might 
insert  magnitude  and  phase  errors  must  be  accounted  for  in  array  design  [11].  The  effects 
on  the  pattern  can  include  loss  in  gain,  increased  sidelobe  levels,  beam  pointing  errors 
and  increased  cross  polarization. 

Many  factors  contribute  to  deviations  in  phase.  The  passive  microwave  circuitry 
in  the  antenna  is  one  source  and  can  be  considered  to  be  independent  of  time.  Other  phase 
errors  (as  well  as  amplitude  errors)  are  introduced  in  the  mixing  process.  Phase  noise, 
phase  offset  and  inaccurate  frequency  are  some  problems  found  in  non-ideal  devices. 
Likewise,  magnitude  error  can  occur  in  any  stage  of  transmission  and  reception  due  to 
gain  mismatch,  interference,  or  noise.  This  type  of  error  can  affect  the  sidelobe  level 
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when  using  taper  to  lower  sidelobes.  Furthermore,  array  elements  are  likely  to  be 
physically  displaced  by  a  fraction  of  wavelength,  which  is  caused  by  array  surface 
deflection,  installation  errors  and  surface  vibration;  thus,  when  elements  are  displaced 
from  their  designed  location,  phase  errors  increase. 


a.  Random  Phase  Noise 


The  phase  noise  represents  frequency  fluctuations,  which  have  their  origin  in 
either  internal  devices  or  external  phenomena,  such  as  mechanical  vibrations  and 
electromagnetic  interference.  For  timing  and  metrology,  instead  of  phase  noise  in  the 
frequency  domain,  one  can  measure  jitter,  which  is  a  time  fluctuation.  For  RF  design 
purposes,  phase  noise  is  preferred  because  it  is  obtained  directly  from  a  spectrum 
analyzer  [15]. 


Phase  noise  is  modeled  as  frequency-selective  random  process  and  is 
characterized  using  the  power-law  model,  which  is  composed  by  five  types  of  noise, 
discussed  further  in  Section  IV.A.3: 

1.  Random  walk  frequency  modulation  noise 

2.  Flicker  frequency  modulation  noise 

3.  White  frequency  modulation  noise 

4.  Flicker  phase  noise 

5.  White  phase  noise 


By  integrating  the  phase  noise  over  the  signal  bandwidth,  one  obtains  the  RMS 
value  for  random  phase  deviation  (p(t) ,  where  (p  is  the  symbol  used  for  phase  noise, 
whereas  4>  is  the  spherical  azimuth  variable,  and 


where 


<Prms  ~ 


-il/2 


J  Sv(f)df 
Jr 


(3.7) 


fa  =  starting  reference  frequency  far  from  signal  of  interest  {fp»  f0),  and 
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phase  noise  spectrum. 


& 


<P 


Notice  the  expression  for  ’s  bandwidth  dependent;  hence,  the  use  of  a 

narrowband  source  is  preferable.  For  example,  the  phase  noise  inserted  by  a  crystal 
oscillator  is  usually  small  because  this  type  of  oscillator  usually  has  a  high  Q  factor  [10]. 

Although  the  random  phase  deviation  <p(t)  is  not  constant  over  time,  if  <p(t)  can 
be  modeled  as  a  wide  sense  stationary  process,  then  at  time  ta,  the  phase  deviation  can  be 
expressed  as  (p{tn)  =  (p,  which  is  independent  from  the  time  origin.  In  this  case  the 

random  phase  model  is  accounted  in  the  array  factor  expression  (3.5)  by  adding  eJ<^>n 5 
which  is  the  phase  error  due  to  phase  noise  at  the  nth  element.  The  modified  expression  of 
array  factor  including  phase  noise  is 


AT 


(3.8) 


n- 1 


b.  Scan  Weight  Errors 


In  (3.5),  the  term  Wn  represents  the  ntb  element  complex  weight.  This  coefficient 
is  precisely  selected  in  order  to  coherently  phase  the  element  transmissions  and 
receptions  to  obtain  the  maximum  gain  for  a  chosen  direction  \f/Sn .  The  magnitude  of  Wn 

is  also  used  to  control  sidelobe  levels  through  tapering.  In  a  digital  array,  the  coefficient 
is  obtained  through  quadrature  modulation. 

The  complex  weight  is  a  part  of  the  array  factor  expression  and  is  expressed  in 
polar  or  rectangular  form  as 


(3.9) 


where  +  ysn.  The  weight  is  upconverted  to  RF  frequency  ^a>Q  =  2 to  the 

in-phase  and  quadrature  signals  given,  respectively,  by 
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(3.10) 


In(t)  =  Rej^e'1'”  jcos(cy )  =  An  cos(Tw)cos(ft>0/) 
and 

Qn(0  =  \m\^AneVn  Jcos(cy)  =  Jwsm(lPM)sm(®0/).  (3.11) 

The  result  after  subtracting  (3.10)  from  (3.11)  is 


S„(t)=  Ancos(a)o(  +  'Y  n). 


(3.12) 


The  expression  for  Sn(t )  represents  the  transmitted  RF  signal  by  the  nth  element. 
The  polar  form  of  the  complex  weight  is  obtained  from  the  relationships 


and 


4,  =  ^.2(O+02(O 

(3.13) 

^  o 

(3.14) 

For  demodulation,  the  received  signal  is  downconverted  to  I  and  Q  channels  using 
a  process  similar  to  transmission.  Let  Rn(t)  be  the  received  RF  signal  by  the  nth  element, 
which  is  expressed  as 


R„(0  =  Kn(t)cos[<o0t  +  *■„(/)]  =  Re |  Kn{i)e 


JKn(  0 


J<°J 1 


(3.15) 


fC  f/) 

where  Kn(t)e  "v  ’  is  the  received  signal  complex  envelope.  A  low  pass  filter  is  applied 

to  reject  the  high  frequency  component  [20].  For  direct  downconversion,  the  expression 
for  In(t)  is 
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I„(t)  =  LPF{Rn(t)cos(a,0t)} 


=  |iw|Re 


f 


Kn(  0 


e/*c«(o+ejK(o+2^) 


=  ^„(')cos[«c„(0] 

where  denotes  a  low-pass  filter  operator.  Similarly, 


Qn(l)  =  LPF{Rn(t)sm(co0t)} 


-LPFlRe\Kn(t) 


ej*n(  0+eiK(0+2®ot) 


(3.16) 


(3.17) 


=  ^„(0sm[AfM(0]- 


The  /„(r)  and  <2«(0  are  sampled  and  recombined  in  the  digital  domain  and  the 
complex  envelope  is  recovered.  Lastly,  the  In(t)  and  Qn(t)  are  multiplied  by  the  complex 
weight  Wn  in  order  to  obtain  the  maximum  radiated  power  in  the  scanned  direction  set  by 
phase  shift  \ps  - 

Any  unknown  magnitude  and  phase  errors  that  modify  In(t)  and  Qn(t )  during  the 
demodulation  process  are  not  incorporated  into  the  complex  weight  coefficient. 
Consequently,  these  errors  change  the  array  factor;  therefore,  one  quadrature  modulation 
and  demodulation  weakness  is  the  required  balance  between  I  and  Q  channels. 

The  following  sources  can  lead  to  significant  phase  and  amplitude  errors  [20]: 

1 .  Phase  error  between  the  I  and  Q  channels 

2.  Differential  gain  error  between  the  two  channels 

3.  Leakage  of  the  LO  signal 

4.  DC  offsets  in  the  I  and  Q  signals 
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An  expression  for  complex  weight  accounting  for  quadrature  errors  is 


^=o+^)V<T',+'") 

where  Sn  is  the  magnitude  error  and  en  is  the  phase  error. 


(3.18) 


Introducing  the  modified  complex  scan  weight  into  array  pattern  expression  of 
(3.5),  we  get 


N 


N 


AF(6,(f>)  =  X  Wnejk'r"  =  X  (\  +  8n)Aneji'V"+e»)ejk'rn . 
11—  1  11— \ 


(3.19) 


c.  Element  Location  Errors 

The  array  factor  term  e^k'r”  is  the  phase  shift  to  element  n.  By  knowing  this 
information,  we  can  set  the  scan  weight  to  cancel  this  phase  effect  in  the  direction  of  the 
beam;  however,  in  case  of  imprecise  location  of  the  element,  this  can  cause  a  phase  error 
that  affects  the  beam  pattern.  If  the  element  position  is  expressed  as  [12] 

rn  =  r0n  +  Ar ’n  =  [(*„  +  Ax)x  +  (yn  +  A y)y  +  (zn  +  A z)z] ,  (3.20) 

the  array  factor  described  in  (3.5)  can  be  rewritten  as 

N  r  _  N  r  -  r  - 
AF(6,4>)  =  X  WnejHr°n+^n)  Wnejk'r°nejk'A,'n .  (3.21) 

11 -\  11— \ 


Note  that  the  phase  error  due  to  element  displacement  is  angle  dependent. 
d.  Final  Expression  for  Array  Factor 

In  summary,  the  phased  array  is  sensitive  to  phase  and  magnitude  errors  caused 
by  random  phase  noise,  quadrature  modulation  and  demodulation  errors,  and  element 
dislocation.  Considering  all  the  possible  phase  and  magnitude  errors,  we  get  the  final 
expression  for  array  factor  as 


30 


(3.22) 


N 


AF(d,$)=J£ 

n- 1 


0  +  <5> 


j(<Pn+en+k'Ar") 


4  e^ne^'r° 


where  yPn+en  +  k»Arnj  is  the  total  phase  error  due  to  phase  noise,  I  and  Q  unbalance 
and  localization  errors,  whereas  (l  +  <^„)  is  magnitude  error. 


5.  Directivity  and  Gain 

The  directivity  is  an  important  characteristic  of  antenna  performance  and  can  be 
expressed  in  terms  of  normalized  pattern  factor  Fnorm(@,$)[3]  : 

An  An 


D{Oj)  = 


Inin 


o  0  \Ft 


max\ 


'In'ln 

u 

0  0 


(3.23) 


J  J  ^  sin 0d0d<f>  J  {  \Fnorn\  sm0d0d^ 


In  the  case  of  a  lossy  array,  the  gain  can  be  computed  from  the  directivity  by  including 
the  antenna  efficiency  T]  [3] : 


G(e,  f) = nixm = —f  [FnrJ.e,^ 

A, 


(3.24) 


where 


T]  =  the  antenna  efficiency,  0  <  r]  <  1, 


PradjPm 

Prad  =  power  radiated  by  the  antenna, 

Pin  =  power  into  the  antenna,  and 
Ae  =  the  array  effective  area. 

The  peak  gain  G0  is  the  maximum  value  of  the  gain  [3],  G0  =  max { (7(#,^)} .  Except  for 
very  simple  antennas,  the  integral  in  (3.23)  must  be  evaluated  numerically. 
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6.  Array  Average  Power  Pattern  and  Sidelobe  Level 

There  are  two  approaches  to  evaluate  the  effect  of  amplitude  and  phase  errors  on 
array  performance.  The  errors  can  be  modeled  as  random  and  Monte  Carlo  simulations 
can  be  performed.  A  large  number  of  independent  trials  are  conducted  and  the  results 
averaged  to  estimate  the  expected  gain,  sidelobe  level  and  beam  pointing  error.  This  is 
the  approach  taken  in  Chapter  V. 

The  array  average  power  pattern  is  equivalent  to  averaging  the  square  of  the 
absolute  value  of  the  array  pattern,  given  by  (3.1),  which  represents  the  radiated  electric 
field.  For  an  array  with  N  elements,  the  average  power  pattern  was  derived  in  [3], 
Assuming  small,  zero  mean  Gaussian  errors,  we  get 

=  (3.25) 

where  {•}  denotes  expected  value.  The  first  term  is  the  error  free  pattern  reduced  in 
intensity.  The  second  term  is  the  average  sidelobe  level. 

The  expected  power  of  the  main  lobe  is 


and  the  average  side  lobe  ratio 


(3.26) 


^sidelobe  j  _  ^ 
^ mainlobe  } 


B.  SUMMARY 


(3.27) 


In  this  chapter  the  digital  array  theory  was  reviewed  by  revisiting  the  expressions 
for  calculation  of  the  array  pattern,  element  factor,  array  factor,  gain  and  sidelobe  level, 
and  modified  equations  were  proposed  for  the  array  pattern,  including  phase  errors, 
which  are  used  in  Chapter  V  simulations.  In  the  next  chapter  the  mathematical  model  for 
instantaneous  phase  for  non-ideal  LO  and  the  random  phase  noise  model,  which  is  used 
in  array  pattern  simulations,  is  presented. 
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IV.  TIMING  AND  SYNCHRONIZATION 


Operation  of  a  distributed  digital  array  relies  on  phase  and  time  synchronization 
of  large  number  of  elements  in  order  to  coherently  process  transmitted  pulses  and  achieve 
high  efficiency.  Even  though  time  and  phase  synchronization  are  closely  related,  these 
concepts  are  applicable  to  different  situations. 

In  order  to  describe  time  synchronization,  first  it  is  necessary  to  define 
chronosignal  or  timing  signal,  which  is  a  signal  used  to  control  actions  that  are  required 
to  be  timed  and  is  generated  by  a  clock  [15].  Ideally,  clocks  with  identical  frequency 
generate  chronosignals  that  differ  only  by  a  constant  offset,  but  in  reality  not  even  the 
world’s  most  important  time  standard  is  perfectly  synchronized  with  the  Earth’s  rotation 
clock  standard.  The  coordinated  universal  time  (UTC)  is  periodically  adjusted  using  leap 
seconds  to  keep  synchronized  with  legacy  standard.  The  leap  second  information  is 
transmitted  for  applications  that  use  UTC  as  a  standard  [21];  therefore,  the  notion  of  time 
synchronization  depends  on  the  application.  The  definition  from  [22]  considers  two 
clocks  to  be  synchronized  if  both  have  the  same  time  origin  and  time  interval  and  do  not 
differ  by  a  specified  uncertainty. 

Phase  is  also  a  linear  function  of  time,  but  phase  has  a  physical  meaning  for 
periodic  signals.  To  perform  a  frequency  demodulation,  the  received  signal  and  local 
oscillator  must  be  phase  synchronized,  i.e.,  locked;  even  if  the  signal  phases  are  360 
degrees  apart,  they  remain  synchronized  because  of  the  periodicity. 

Every  distributed  sensor  network  requires  time  synchronization,  and  many 
improvements  have  been  achieved  in  recent  years  to  provide  better  performance. 
Likewise,  a  digital  array  radar  is  a  type  of  sensor  network,  consisting  of  nodes  that 
transmit  and  process  radar  pulses  coherently.  A  digital  array  also  demands  phase 
synchronization  for  beamforming  purposes. 
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A. 


TIME  AND  PHASE  SYNCHRONIZATION 


A  mathematical  model  for  timing  signals,  generally  used  for  clock  and  oscillator 
characterization,  is  described  in  this  section.  The  type  of  oscillator  depends  upon  the 
design  and  application.  There  are  a  variety  of  devices  available  on  the  market  using 
different  technologies  such  as  those  described  in  the  Appendix. 

In  [15]  the  author  divides  clocks  in  two  groups:  autonomous  and  slave  clocks. 
Autonomous  clocks  are  those  that  have  an  internal  oscillator  and  do  not  depend  on 
external  influence  to  supply  a  timing  signal,  such  as  crystal  oscillators;  slave  clocks  are 
those  that  require  external  control. 

1.  Frequency,  Phase  and  Time  Models 

The  model  used  in  [15]  to  depict  the  behavior  of  instantaneous  frequency  in  real 
clocks  consists  of  three  terms  given  by 

m=fs+M')+fA<)  (4.D 


where 

fs  =  starting  frequency, 

fn  +  /offset 

fn  =  nominal  frequency, 

/  t  =  frequency  offset, 

fd(t)  =  frequency  drift,  and 
fa(t)  =  random  frequency. 

This  definition  is  not  specific  for  a  particular  oscillator  or  signal.  Any  synthesized 
continuous  wave  (CW)  signal  is  subject  to  these  phenomena.  Although  the  focus  here  is 
the  local  oscillator  frequency  errors  fLO,  (4.1)  is  also  applied  for  errors  in/0  and  f0iF. 
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The  starting  frequency  is  a  constant  value  and  includes  the  oscillator  nominal 
frequency  plus  an  initial  calibration  error;  in  other  words,  it  is  the  device  initial  accuracy. 
Another  factor  that  affects  the  instantaneous  frequency  is  the  drift,  which  is  a 
time-dependent  factor,  and  older  devices  tend  to  be  less  accurate.  Usually,  frequency  drift 
is  approximated  by  a  linear  function,  called  linear  fractional  frequency  drift  rate  [15]: 

m=Drj  (4,2) 

where  D  is  the  linear  fractional  frequency  drift  rate. 

The  last  term  in  (4.1)  is  a  random  frequency  deviation  that  represents  the 
oscillator  phase  noise 


fa(,)  2jt  1  <43) 

where  the  term  (p  represents  the  random  phase  deviation,  and  the  instantaneous 

frequency  expression  can  be  rewritten  as 


/<o =/,+%*+ 


1  dq> 
2  n  dt 


(4.4) 


Consequently,  from  (4.4)  it  is  possible  to  derive  an  expression  for  instantaneous 


phase 


+  +(p(i)-\-Q>0  (4.5) 

where  (D0  is  the  initial  phase.  Finally,  the  expression  for  a  time  signal  generated  by  an 
oscillator  is 


7X/)  = 


2  */„' 


Note  that  for  the  ideal  case,  where  D  =  0 ,  <p(t)  =  0  and  <D0  =  0, 


7X0  = 


m 

27TVn 


2  7TV„t 
—  =t 
2jrvt j 


(4.6) 


(4.7) 
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2.  Accuracy  and  Drift 

The  maximum  frequency  error  of  a  clock  during  its  life  cycle  determines  the 
accuracy.  Frequency  error  is  mainly  affected  by  the  offset  and  drift  and  is  expressed  as  a 
ratio  of  clock  deviation  to  its  nominal  frequency  [15] 

f€  =  Mmi  (4.8) 

Jn 

where  Afmax  is  the  maximum  frequency  deviation. 

The  oscillator  frequency  drifts  due  to  internal  or  external  phenomena.  The  main 
internal  factor  is  the  aging  of  the  oscillator  components;  although,  it  is  very  difficult  to 
determine  an  exact  reason  for  internal  drift  [15]. 

Drift  can  also  occur  by  external  factors  such  as  interference,  temperature  and 
vibration.  Oscillators  in  modulators  and  demodulators  are  subject  to  injection  pulling  and 
locking,  i.e.,  a  frequency  disturbance  due  to  influence  of  an  external  signal  which  is  close 
in  frequency  to  that  second  oscillator  or  device.  When  the  oscillator  starts  to  resonate 
with  the  frequency  of  the  disturbance,  this  is  called  injection  locking,  but  if  two 
frequencies  remain  different,  this  is  called  injection  pulling  [23], 

Frequency  accuracy  is  essential  for  distributed  array  functioning.  A  slight 
frequency  deviation,  after  a  long  period  of  time,  might  represent  a  significant  phase 
difference  that  degrades  pattern  characteristics. 

3.  Phase  Noise 

Phase  noise  is  composed  by  the  oscillator  frequency  and  amplitude  fluctuations 
affecting  other  components  that  use  this  signal  as  reference.  There  are  other  sources  of 
phase  noise,  but  oscillators  tend  to  increase  the  noise  close  to  its  resonance  and 
harmonics  [24],  For  this  reason,  after  the  mixing  process,  oscillator  noise  is  more 
significant  than  other  contributors. 

Phase  noise  can  be  characterized  in  the  time  domain  or  the  frequency  domain. 
Time-domain  models  observe  deviations  during  a  period  of  time  and  establish  a  statistical 
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model  by  measuring  the  Allan  variance,  for  instance.  In  metrology,  time-variance 
measurements,  such  as  the  Allan  variance,  are  preferable  because  long-term  variations 
are  registered.  Furthermore,  frequency-domain  model  parameters  are  estimated  from 
time-domain  observations,  which  also  include  drift  effects  [15].  The  model  used  in  this 
thesis  assumes  independent  phase  noise  and  drift. 

Frequency-domain  models  have  been  used  for  decades  and  accurately  represent 
oscillator  noise  over  short  time  intervals.  Moreover,  frequency-domain  phase  noise  can 
be  obtained  directly  from  a  spectrum  analyzer,  simplifying  oscillator  design  and 
characterization  [10]. 


a.  Power-Law  Model 


To  characterize  the  oscillator  phase  noise,  the  power-law  model  is  usually  used, 
which  is  expressed  as  one-sided,  power  spectral  density  (PSD)  of  random  time  deviation 
r(l)  [15] 


ST(f)  = 


1 


o 


2  £  *a+2/“.  0 <f<fh 

a—— 4 

f>fh 


(271)  a=_4 

0, 


(4.9) 


where  /l_2,  h  h$,h\,  hj  are  device-dependent  coefficients.  Each  particular  clock  is 
defined  by  a  set  of  coefficients.  The  term  is  the  oscillator  bandwidth  cut-off 
frequency.  The  relationship  between  random  time  deviation  and  random  phase  deviation 
is 


r(/)  = 


2  xfn 


Consequently,  the  power  spectra  are  related  by 


ST(f)  = 


S9{f) 

(2  nfnf 


(4.10) 


(4.11) 


37 


The  phase  noise  power  is  related  to  the  RMS  value  of  the  random  phase  deviation 
*PRMS ,  expressed  in  (3.7).  The  parameters  h. 2,  h.  1,  h0,  hi,  are  related  to  five  different 
types  of  noise  components: 

1.  Random  walk  frequency  modulation  [  h_2  //4  ] 

Random  walk  frequency  modulation  (RWFM)  is  the  part  of  phase  noise  attributed 
to  environmental  perturbation  such  as  mechanical  shocks,  temperature  variation  [15],  and 
random  changes  in  oscillator  frequency  [25].  RWFM  has  main  frequency  components 
close  to  oscillator  resonance.  For  this  reason  it  is  more  significant  near  these  frequencies 
and  also  more  difficult  to  measure  [15]. 


2.  Flicker  frequency  modulation  |^/?j  //3  J 


The  source  for  flicker  frequency  modulation  (FFM)  is  attributed  to  a  physical 
resonance  mechanism  of  an  oscillator  or  issues  controlling  electronic  components.  Some 
oscillators  might  not  show  FFM  noise,  which  might  be  covered  by  white  frequency 
modulation  (WFM)  or  FPM,  but  in  high  quality  factor  oscillators,  FFM  is  more 
significant. 


3.  White  frequency  modulation  ^hQ  / f 2  J 


White  frequency  modulation  is  generated  by  electronic  components,  and  for  most 
oscillators,  white  frequency  noise  is  less  significant  close  to  the  resonance 
frequency  [25], 

4.  Flicker  phase  modulation  [v/] 

Flicker  phase  modulation  (FPM)  is  a  colored  type  of  noise  and  is  related  to 
oscillator  resonance  mechanism.  In  mixing  and  amplification,  the  local  oscillator  noise  is 
amplified,  especially  close  to  carrier  frequencies  [15]. 

5.  White  phase  modulation  [*2] 

White  phase  modulation  (WPM)  noise  is  a  wideband,  uncorrelated  noise  with 
constant  magnitude  over  the  spectrum  and  can  be  interpreted  as  a  phase  noise  floor.  The 
main  source  for  WPM  is  the  noise  inserted  by  electronic  devices.  Even  though  oscillators 
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have  low  WPM  noise,  circuits  that  employ  digital  devices  increase  WPM  due  to  errors  in 
clock  reference  and  phase  quantization  errors  [15]. 

The  PSD  of  random  time  deviation  in  its  components  is  shown  in  Figure  13,  and 
an  example  of  phase  noise  from  VCO  manufacturer  Minicircuits  is  shown  in  Figure  14. 


Figure  13.  Random  Time  Deviation  PSD  Obeying  the  Power-Law  Model 

Plotted  on  a  Log-Log  Diagram 

The  frequency  /R  is  a  reference  frequency  far  from  the  carrier  frequency.  The  figure  of 
interest  is  the  area  under  the  phase  noise  curve,  which  represents  the  power  and  is  related 
to  random  phase  deviation.  Adapted  from  [15]:  S.  Bregni,  Synchronization  of  Digital 
Telecommunications  Networks,  Vol.  1.  John  Wiley  &  Sons,  2002,  p.  248. 
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ROS-2600-1 1 19+ 
Phase  Noise 


1  10  100  1000 
Frequency  Offset  (KHz) 

Figure  14.  Phase  Noise  of  the  Voltage  Controlled  Oscillator  Model  ROS- 

2600-1119,  Manufacture  Mini  circuits 

Source  [26]:  (2015,  Nov.  10)  Voltage  Controlled  Oscillator  ROS-2600-1119+.  [Online]. 
Available:  http://www.minicircuits.com/pdfs/ROS-2600-l  1 19+.pdf. 


b.  Leeson’s  Equation 

In  1966,  Leeson  [27]  showed  “a  heuristic  derivation,  presented  without  formal 
proof’  of  the  phase  noise  for  an  oscillator.  The  linear  time-invariant  theory  assumes  the 
oscillator  is  a  linear  device  with  some  correction  to  account  for  its  non-linearities.  His 
first  expression  for  phase  PSD  is  [28] 

„  f  nfkbt, ) 

■VI  p  I 

X  ‘out  } 

where 

Np  =  an  empirical  noise  factor, 

Q  =  oscillator  open-loop  unloaded  Q-factor, 

fQ  =  carrier  frequency  (RF), 


,+(4)J 


(4.12) 
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fR  =  reference  frequency  in  signal  bandwidth, 

Pout  =  output  power, 

Kb  =  Boltzmann's  constant  (  =  1 .38  x  10  23  JK  ’),  and 

Ts  =  system  temperature  (in  Kelvin). 

In  reference  [28]  the  author  describes  a  modified  model  including  the  flicker  noise 
effects.  The  phase  noise  is  usually  given  as  a  dB  quantity 


A/m)  =  10  log 


dBc/Hz. 


(4.13) 


The  simplicity  of  this  model  is  the  great  advantage  over  the  power-law.  While  the 
power-law  model  requires  measurements  to  estimate  the  coefficients  h  2,  h  y,  /zq,  hy,  hj , 
Leeson’s  equation  accounts  for  other  oscillator  parameters  that  are  easier  to  obtain; 
moreover,  this  model  includes  frequency  dependence  for  the  phase  noise.  In 
upconversion,  the  ratio  faj fm  increases  the  noise,  adding  a  factor  of  20  dB/decade, 
which  is  significant  for  high  frequency  radar  applications. 


For  example,  in  [29]  the  author  gave  an  example  of  a  LO  signal  in  the  range  of 
5  GHz  which  was  generated  from  a  10-MHz  crystal  oscillator  using  a  frequency 
multiplication.  This  process  increases  the  phase  noise  by 


£( 5  GHz)  -  (20  dB/decade  K 1 0  log 500 )  +  £(!()  MHz) 

=  (20  dB/decade)26 .99  +  £(10  MHz)  (4. 14) 

=  53.9  +  £(10  MHz)  dB. 

Consequently,  the  RMS  phase  noise  increases  by  245,470. 
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B.  TIME  AND  PHASE  SYNCHRONIZATION  TECHNIQUES 

Applications  that  require  coherent  detection  and  transmission,  like  phased  array 
antennas,  need  carrier,  bit  and  word  synchronization  [30].  In  phase-coded  CW  radars, 
these  three  levels  can  be  translated  as  carrier,  subcode  and  code  synchronization. 

The  most  critical  is  the  carrier  synchronism  because  it  requires  errors  smaller  than 
a  fraction  of  cycle,  which  can  be  sub-nanosecond  for  microwaves.  Code  and  subcodes  are 
more  resilient  to  synchronization  errors  than  the  carrier.  For  example,  a  subcode  period  of 
20  ns  is  used  to  achieve  a  3.0-m  range  resolution,  whereas  a  code  period  is  260  ns  if  a 
Barker  code  modulation  of  length  13  is  used.  Nevertheless,  all  three  levels  of 
synchronization  are  important  to  coherent  transmission  of  the  pulse  and  further  detection. 

In  order  to  limit  our  discussion,  time  synchronization  as  a  high  level  requirement 
is  considered  to  assure  code  and  subcode  coherence;  whereas,  phase  synchronization  is 
required  for  a  phased  array  beamforming.  The  time  estimate  can  be  extracted  through 
some  communications  protocols  such  as  network  time  protocol  (NTP)  or  precision  time 
protocol  (PTP). 

1.  Phase  and  Frequency  Synchronization  Using  Phased  Locked  Loop 

For  a  distributed  array,  two-step  synchronization  is  proposed:  frequency  and 
phase  synchronization.  In  communications,  coherent  detection  also  requires  carrier 
frequency  recovery  either  from  a  pilot  channel  or  through  received  signal  [30],  The 
former  is  equivalent  to  a  direct  distribution  LO  introduced  in  Chapter  II.B.l,  and  the 
latter  uses  a  VCO.  For  carrier  extraction  from  a  received  signal,  a  PLL  is  used  to  control 
the  demodulator  local  oscillator;  this  architecture  is  known  as  a  Costas  loop.  Once  the 
local  oscillator  is  synchronized  to  the  received  signal,  coherent  quadrature  demodulation 
can  be  performed.  We  briefly  introduce  the  important  PLL  concepts  and  operating 
modes. 


a.  Phased  Locked  Loop 

A  PLL  is  a  type  of  slave  clock  that  generates  a  timing  signal  controlled  by  a 
reference  input.  PLLs  are  widely  employed  in  frequency  synchronization.  Despite  the 
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name,  PLLs  are  also  used  as  non-linear  filters  and  in  modulation,  demodulation  and 
frequency  synthesis  [15].  The  basic  loop,  illustrated  in  Figure  15,  is  composed  of  the 
following: 

1.  Phase  Detector 

2.  Loop  Filter  /  Amplifier 

3.  VCO 


Figure  15.  Block  Diagram  of  PLL 

Adapted  from  [15]:  S.  Bregni,  Synchronization  of  Digital  Telecommunications  Networks, 

vol.  1 .  John  Wiley  &  Sons,  Ltd,  2002. 

The  phase  detector  (PD)  is  a  non-linear  circuit  that  compares  the  reference  signal 
phase  and  the  output  from  the  VCO.  If  the  VCO  frequency  is  close  to  the  reference  input, 
then  a  low  frequency  output  is  expected  from  the  PD.  A  low  pass  filter  is  used  to 
eliminate  PD  noise  and  harmonics,  and  its  output  is  feedback  to  the  VCO.  When  the 
frequency  is  locked,  the  VCO  input  becomes  a  DC  voltage. 

A  PLL  has  three  operation  modes: 

1.  Locked  Mode:  In  locked  mode  the  PLL  is  tracking  the  reference  frequency.  For 
an  ideal  case,  the  locked  mode  is  receives  a  stable  and  continuous  reference  signal,  but  in 
real  applications,  noise,  interference  and  component  malfunctioning  can  corrupt  the  PLL 
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input.  For  this  reason  the  locked  mode  operates  in  stressed  mode,  which  means  the  PLL 
“incorporates  phase  build-out  features  to  limit  time-keeping  error  on  transients  after 
micro-interruptions”  [15]. 

2.  Free-run  mode:  In  free-run  mode,  the  PLL  operates  as  an  oscillator,  supplying 
an  output  signal  and  disregarding  the  input  reference.  In  fact,  it  is  very  rare  that  a  PLL 
operates  in  free-run  mode;  it  is  preferable  to  run  it  in  hold-over  mode. 

3.  Hold-over  mode:Similar  to  free-run,  when  the  PLL  is  in  hold-over  mode,  it 
generates  a  signal  autonomously  when  the  input  reference  fails;  however,  a  PLL  in  hold¬ 
over  mode  keeps  constant  the  voltage  that  feeds  the  VCO,  using  the  last  value  before  the 
reference  interruption.  A  PLL  might  even  store  several  past  values,  and  in  case  of 
interruption  the  VCO  is  controlled  using  some  estimation  of  these  values. 

b.  Costas  Loop 

The  Costas  loop  is  a  PLL  application  for  carrier  frequency  and  phase  recovery. 
The  block  diagram  is  illustrated  in  Figure  16.  The  loop  is  a  coherent  quadrature 
demodulator  with  a  VCO.  The  operation  is  similar  to  a  PLL;  however,  the  difference  is 
that  the  output  signal  is  not  the  oscillator  reference,  but  the  demodulated  signal. 


Figure  16.  Block  Diagram  of  Demodulator  using  Costas  Loop 

Adapted  from  [30]  B.  P.  Lathi,  Modern  Digital  and  Analog  Communication  systems ,  2nd 
ed.  Sacramento,  CA:  The  Dryden  Press,  1989. 
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c.  Synchronization  Using  PLLs 

The  DBF  transmits  the  control  signals  to  elements  through  a  RF  channel,  which 
can  be  wired  or  wireless.  A  PLL  can  recover  the  carrier  frequency  and  use  it  as  a 
reference  to  generate  a  LO  signal,  which  might  have  a  frequency  higher  than  its 
reference.  The  reason  is  because  radars  usually  operate  at  frequencies  in  the  S,  X  and  K 
bands,  whereas  the  WLAN  might  be  at  2.4  GHz. 

A  Costas  loop  at  each  element  recovers  the  frequency  and  phase  of  the  received 
signal;  however,  the  carrier  cannot  be  used  directly  as  LO  because  the  link  distance  adds 
a  time  delay,  which  cannot  be  neglected.  The  elements  are  subjected  to  different  delays 
because  of  their  various  distances  to  the  DBF.  The  use  of  phase  shifting  to  compensate 
the  delay  was  experimented  by  Loke  [12]  and  Tsai  [13]  in  a  brute  force  synchronization 
approach,  which  achieved  synchronization  phase  errors  as  low  as  22.5  degrees. 

In  [12]  and  [13],  the  authors  used  a  circulator  to  isolate  the  received  broadcasted 
LO  signal  from  the  element  LO  signal  that  is  returned  to  the  DBF  for  comparison  with  a 
master.  The  PLL  synchronization  eliminates  this  device  since  the  LO  is  extracted  from 
WLAN  signal,  and  an  independent  path  is  used  to  transmit  back  the  LO  signal  to  the 
DBF. 

Once  the  phase  shifters  are  set  and  the  PLL  is  locked  to  the  reference,  the  element 
synchronization  is  accomplished.  Theoretically,  the  element  might  not  lose 
synchronization  unless  the  frequency  of  the  DBF  reference  carrier  changes. 

There  is  a  limitation  with  this  technique.  The  WLAN  or  wired  LAN  must  have  a 
single  carrier  frequency.  Frequency-hopping  channels,  like  WLAN  standards  IEEE 
802. 1 1  n/g,  might  not  provide  the  PLL  a  stable  LO. 

2.  Phase  and  Frequency  Synchronization  Stable  Local  Oscillator 

A  STALO  can  be  used  for  phase  and  frequency  synchronization  of  a  digital 
distributed  array.  One  example  of  a  STALO  is  a  chip  scale  atomic  clock  (CSAC)  [31], 
which  is  a  single  frequency  clock  reference  that  can  be  used  directly  with  a  digital  signal 
processing  device,  like  a  DDS,  or  feeding  a  PLL  to  synthesize  a  high  frequency. 
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With  the  CSAC,  phase  synchronization  can  be  achieved  using  digital  delays  once 
the  output  signal  is  a  square  shape  waveform.  The  CSAC  also  provides  a  pulse  per 
second  (PPS)  input  for  synchronization  which  can  be  used  after  compensating  for  the  link 
distance  delay. 

C.  SUMMARY 

In  this  chapter  a  mathematical  model  for  non-ideal  LOs  which  have  phase  noise 
and  inaccurate  frequency  was  presented.  Additionally,  LO  distribution  techniques  that 
can  minimize  phase  errors  were  presented. 

Results  from  parametric  simulations  programmed  in  MATLAB  are  presented  in 
Chapter  V,  with  the  objective  of  illustrating  the  effects  of  phase  errors  on  pattern 
distortions  and  to  aid  the  array  design  process  to  set  acceptable  limits  for  these  types  of 
errors.  These  simulations  were  based  on  the  models  and  equations  presented  in  this 
chapter  and  Chapter  III. 
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V.  SIMULATIONS  AND  RESULTS 


A.  INTRODUCTION 

A  digital  array  beamformer  controls  each  array  element  to  generate  the  waveform 
and  radiation  pattern.  The  fact  that  all  elements  operate  individually  allows  for  the 
correction  of  phase  and  magnitude  errors  that  originate  during  implementation  [3]. 

There  are  random  errors  that  are  very  difficult  to  eliminate  or  compensate.  In 
large  arrays  these  errors  can  limit  the  minimum  sidelobe  level  and  maximum  array  gain 
[32],  The  block  diagram  in  Figure  17  depicts  the  T/R  module  proposed  in  [3]  for  a 
distributed  array,  and  it  has  highlighted  in  red  the  main  sources  of  phase  and  amplitude 
errors  that  affect  the  array  pattern  which  were  considered  in  this  study:  LO  phase  noise, 
LO  frequency  linear  drift  and  I  and  Q  gain  unbalance. 


Figure  17.  A  T/R  Module  Block  Diagram.  Phase  Error  Sources  are 

Highlighted  in  Red. 

Adapted  from  [3]:  D.  Jenn,  Y.  Loke,  T.  C.  H.  Matthew,  Y.  E.  Choon,  O.  C.  Siang,  and  Y. 
S.  Yam,  “Distributed  Phased  Arrays  and  Wireless  Beamforming  Networks,”  Int.  J. 
Distrib.  Sens.  Networks,  vol.  5,  no.  4,  pp.  283-302,  2009.  p.  284. 
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The  LO  supplies  the  reference  signal  for  the  modulator,  demodulator  and  all 
digital  components  that  demand  a  clock  reference.  Phase  errors  are  generated  if  LOs 
errors  are  propagated  into  these  components.  A  wireless  LO  distribution  concept  is 
shown  in  Figure  17;  however,  most  applications  employ  a  hardwired  distribution  circuit. 

In  Chapter  III.A.4.b,  the  consequences  of  unbalanced  gain  between  the  I  and  Q 
branches,  which  adds  a  phase  error  that  can  degrade  the  array  pattern  characteristics,  was 
explained.  Mislocation  of  array  elements  is  another  source  of  phase  error.  Depending  on 
the  operating  frequency,  a  couple  of  centimeters  can  represent  more  than  one  wavelength. 

These  phase  errors  were  modeled  using  MATLAB,  and  the  simulation  results  are 
presented  in  the  following  sections. 

B.  ARRAY  FACTOR  COMPLEX  WEIGHTS 

The  array  factor  expressed  in  (3.5)  has  the  complex  weight  Wn,  which  is  set  to 
scan  the  main  beam,  control  sidelobes  and  also  potentially  compensate  for  phase  errors 
[13].  In  this  simulation  scenario,  the  magnitude  of  the  complex  weight  is  used  for 
tapering,  whereas  the  phase  is  employed  for  beam  scanning. 

Phased  array  antennas  designed  for  radars  typically  use  a  Taylor  taper  for  si  delobe 
level  control.  The  main  characteristics  of  arrays  that  use  this  taper  are  the  monotonically 
decreasing  sidelobes  and  the  absence  of  edge  discontinuities  [33]. 

For  the  linear  array  simulations,  the  taper  was  set  to  obtain  a  maximum  relative 
sidelobe  level  of  -35  dB  and  Nbar  of  5,  which  is  the  number  of  “nearly  constant-level 
sidelobes  adjacent  to  the  mainlobe”  [34],  For  the  array  of  subarrays  simulations,  for  each 
subarray  the  taper  was  set  to  a  maximum  relative  sidelobe  level  of  -20  dB  and  an  Nbar  of 
2.  For  beam  pointing,  the  scan  angles  are  given  by  direction  cosine  expressions 
us,  vv ,  Wy  in  (3.1),  nulling  the  exponent  in  the  scan  direction. 

C.  PATTERN  CHARACTERISTICS 

The  purpose  of  these  simulations  is  the  numerical  calculation  of  the  array  pattern 
characteristics,  i.e.,  gain,  average  and  peak  sidelobe  levels,  and  pointing  errors.  The  array 
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pattern  formula  given  in  (3.1)  was  programmed  in  MATLAB  to  simulate  hardware  errors 
in  the  arrays,  and  the  results  were  used  to  extract  the  sidelobe  information  and  beam 
pointing  direction. 

1.  Average  sidelobe  level 

The  average  SLL  is  computed  from  pattern  samples  from  outside  the  mainbeam 
as 


{SLL}  =  ~—  f;  (5,1) 

n=-N  n-n2 

where 

(SLL)  -  average  sidelobe  level, 

F{0,<f>)  =  array  pattern, 

SO  =  numerical  calculation  step  size  (sampling  interval)  in  degrees, 

N p  =  Array  maximum  angle,  generally  90°, 


«l  =  mainlobe  first  left  null,  and 

«2  =  mainlobe  first  right  null. 


For  this  simulation  <j>  is  assumed  to  be  zero,  and  the  step  size  \0  is  0.2  degrees.  The 
mainlobe  first  nulls,  left  and  right,  were  found  using  MATLAB  function  findpeaks. m. 

2.  Peak  sidelobe  level 

In  similar  fashion,  the  peak  sidelobe  level  relative  to  the  mainbeam  maximum  is 
expressed  by: 

MAX(SLL)  =  MAX^F(nSO,<f>f^,  for/i  <«l  andtt  >n2  (5.2) 

where  MAX(-)  represents  the  maximum  value  of  the  vector. 

3.  Beam  pointing  angle 

The  beam  pointing  direction  is  the  angle  0max  at  which  the  radiation  pattern  has 
the  maximum  value 
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(5.3) 


FiO^j)  =  MAX^F{n/\(h<t>f^  fern,  <n<  n2. 

Ideally,  the  value  for  6^^.  is  equal  to  0  ,  the  scanning  angle;  therefore,  the  difference 
between  the  0max  and  0  is  the  beam  pointing  error 


BPE  =  0max-0s.  (5.4) 

Under  the  influence  of  random  phase  noise,  random  7  and  Q  gain  unbalance,  and  element 
mislocation,  the  BPE  is  a  random  variable  which  has  a  standard  deviation  relative  to  the 
half-power  beamwidth 


aBPE~ 


(  I VA /?((), 


max 


0: 


HPBW 


(5.5) 


where  the  VAR(-)  is  the  variance  and  the  Ojjpgjy  ’s  the  half-power  beamwidth.  In  this 
thesis,  the  mean  of  the  BPE  is  not  expressed  relative  to  0/jpp^y  but  is  given  by 

( BPE)  =  {9max-9s ).  (5.6) 

4.  Pattern  gain 


The  gain  is  calculated  using  (3.24),  which  requires  the  normalized  antenna 
radiation  pattern  integration;  therefore,  the  Gaussian  quadrature  technique  was  employed. 


Gaussian  quadrature  is  a  numerical  integration  technique  whose  result  approaches 
the  optimum  integration  value  [35],  It  uses  a  weighting  function  and  selects  abscissa 
values  to  match  the  roots  of  a  preselected  orthogonal  polynomial  which  depend  on  the 
integration  interval.  The  simulation  abscissas  and  weighting  function  can  be  found  in 
many  handbooks  of  mathematical  functions,  for  example,  reference  [36], 


D.  UNIFORM  ARRAY  HARDWARE  ERROR  SIMULATIONS 

1.  Simulation  Scenario  for  Phase  Noise,  Gain  Unbalance  and  Frequency 
Drift  Error 

In  Table  1  the  simulation  settings  for  a  uniform  linear  array  of  100  identical 
elements  employed  in  the  phase  noise,  7  and  Q  gain  unbalance  and  frequency  linear  drift 
simulation  scenarios  is  summarized. 
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Table  1.  Array  Characteristics  used  in  MATLAB  Simulations 


Number  of  elements 

100 

Elements  spacing 

0.2  wavelength 

Element  factor 

cos(0) 

Nominal  frequency 

300  MHz 

Tapering 

Taylor  (-35  dB  SLL,  Nbar  =  5) 

2.  Phase  Noise  Simulation 

a.  Simulation  Scenario 

Non-ideal  LOs  are  subject  to  phase  noise,  i.e.,  amplitude  and  phase  fluctuations 
caused  by  internal  and  external  factors.  The  phase  noise  characterization  is  detailed  in 
Chapter  IV. A.  3.  For  simulation  purposes,  the  phase  noise  was  modeled  as  WPM 
converted  to  time-domain  as  random  phase  noise  deviation  q),  which  is  applied  in  the 
modified  expression  for  array  pattern  given  in  (3.8).  The  estimation  of  q)  is  addressed  in 
Section  F.  1 . 

This  simulation  scenario  assumes  independent  LOs  for  all  elements  but  assumes 
they  are  synchronized  with  each  other.  In  other  words,  all  LOs  have  the  same  frequency 
and  phase  offset,  but  the  phase  noise  is  present  and  independent  for  all  oscillators. 

b.  Simulation  Results 

The  simulation  results  are  obtained  from  the  average  of  100  trials.  The  Taylor 
taper  is  set  to  give  a  maximum  relative  sidelobe  of  -35  dB.  Four  different  scanning  angles 
were  simulated:  0°,  15°,  30°  and  45°. 

Shown  in  Figure  18  is  the  pattern  distorted  by  phase  deviation.  The  top  pattern  is 
based  on  a  phase  noise  standard  deviation  of  0.05  rad,  whereas  for  the  bottom  pattern  the 
noise  has  0.5  rad  of  standard  deviation.  The  error  free  pattern  is  shown  by  the  dotted  line. 
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Figure  18.  Linear  Array  Pattern  Degraded  by  Phase  Noise  (solid), 
(a)  <Prms  =  °-05  rad  and  (b)  <Prms  =  0.5  rad 


From  Figure  18,  there  is  a  noticeable  increase  in  sidelobes,  which  is  confirmed  by 
(SLL)  and  MAX  (SLL)  results.  The  plot  in  Figure  19  illustrates  the  average  sidelobe 
level,  calculated  using  (5.1),  as  a  function  of  phase  deviation.  Note  that  for 
VRMS  >0.12  rad,  (SLL)  is  already  larger  than  the  maximum  sidelobe  level  limited  by 
the  Taylor  taper  and  is  worse  for  oblique  scanning  angles. 

In  Figure  20  the  peak  sidelobe  level,  computed  using  (5.2),  as  a  function  of  phase 
deviation  is  plotted.  Note  that  for  (Prms  >0.16  rad,  the  MAX(SLL)  is  ten  times  larger 

than  maximum  desired  sidelobe  level,  which  is  -35  dB.  As  expected,  the  mainlobe  gain 
decreases  as  a  consequence  of  the  sidelobe  level  increase.  In  Figure  21  the  array  gain  loss 
is  illustrated  as  function  of  phase  deviation.  The  loss  is  around  0.5  dB  for 
<PRMS  ~  rad. 
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Figure  19.  Relative  Average  Sidelobe  Level  (SLL'jvs.  Phase  Noise  Standard 
Deviation  <Prms  Plotted  for  Four  Scanning  Angles 


Figure  20.  Relative  Peak  Sidelobe  Level  MAX(SLL)  vs.  Phase  Noise 
Standard  Deviation  (pR^s  Plotted  for  Four  Scanning  Angles 
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Figure  21.  Array  Gain  Loss  Due  to  Phase  Noise  Plotted  for  Four  Different 

Scanning  Angles 


Pointing  error  is  another  concern  for  arrays  performing  tracking.  In  Figure  22  the 
average  beam  pointing  error  in  degrees  is  illustrated,  and  the  beam  pointing  error 
standard  deviation  relative  to  the  HPBW  is  shown  in  Figure  23.  Note  that  for  the  pointing 
error,  standard  deviation  is  larger  than  5%  of  the  HPBW,  which  is  a  significant  error  for 
high  accuracy  systems.  Long-range  applications  such  as  satellite  communications  and 
radar  are  sensitive  to  small  pointing  errors  that  result  in  missing  the  target  by  several 
meters  or  kilometers.  The  HPBW  as  a  function  of  gain  unbalance  is  shown  in  Figure  24. 
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Figure  22.  Average  Beam  Pointing  Error  ( BPE }  vs.  Phase  Noise  Deviation 
(Prms  Plotted  for  Four  Different  Scanning  Angles 


Figure  23.  Relative  Standard  Deviation  <r^pp  vs.  Phase  Noise  Deviation 
(Prms  Plotted  for  Four  Different  Scanning  Angles 
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Figure  24.  Half-Power  Beamwidth  vs.  Phase  Noise  Deviation  (pjurfs  Plotted 

for  Four  Different  Scanning  Angles 


3.  I  and  Q  Gain  Unbalance 

a.  Simulation  Scenario 

This  simulation  scenario  reproduces  the  effects  of  non-ideal  quadrature 
modulators  and  demodulators  which  are  used  in  radar  and  communications  systems  for 
phase  control.  In  Figure  25  a  model  of  the  simulation  scenario  is  illustrated.  Here,  each 
element  is  supposed  to  have  a  random  attenuator  for  7  and  Q  branches,  L/  and  Lq.  These 
attenuators  represent  the  gain  unbalance,  and  their  values  are  normally  distributed.  This 
random  variable  probability  density  function  (PDF)  is  modeled  as 

L2  } 

- y  (5.7) 

W) 

where 

L  =  random  attenuator  loss  (either  L  j  and  Lq ), 

<7  j  =  random  attenuator  standard  deviation. 


PDF  = 


exp 


3/ 2.575,  and 
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d 


nominal  /  and  Q  gain  imbalance. 


Using  the  expression  above  for  crL ,  we  expect  that  the  random  attenuation  will 
be  within  the  nominal  gain  unbalance  values  with  99%  confidence.  The  modified 
expression  for  the  array  factor,  neglecting  the  other  phase  error  sources,  is 


where 

I(Li) 


N  r 


AF(0,$)=2j 
11—  1 


(1  +  8  n  )e/c«  Ae 


,JV 


Jk' 


=i[wnyk-r°n 


11— \ 


^[l(Lj)  +  jQ(LQ)\e’k‘r 

n= 1 


(5.8) 


=  in-phase  signal  distorted  hy  random  attenuation  L /,  and 


»Lq) 


quadrature  signal  distorted  hy  random  attenuator  Lq. 


Figure  25.  Receiving  Simulation  Scenario  where  Lj  and  Lq  are  Hypothetical 

Random  Attenuators 
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b. 


Simulation  Results 


Similar  to  previous  cases,  the  simulation  results  were  obtained  from  the  average 
of  100  trials  for  four  different  scanning  angles:  0°,  15°,  30°  and  45°. 

The  pattern  distorted  by  I  and  Q  gain  unbalance  is  shown  in  Figure  26.  On  the  top 
pattern,  the  elements  have  random  magnitude  attenuation  within  the  limit  of  0.4  dB, 
whereas  for  the  bottom  pattern  the  limit  was  1.0  dB.  The  dotted  line  is  the  pattern  with  no 
errors.  In  Figure  27  and  Figure  28  the  average  sidelobe  level  and  the  peak  sidelobe  level, 
respectively,  are  illustrated.  The  impact  of  I  and  Q  gain  unbalance  is  less  significant  than 
the  phase  noise;  however,  for  L  >  1  dB  this  error  contribution  is  more  noticeable.  A 
L  >  1  dB  is  unlikely  to  happen  because  modulators  and  demodulators  are  designed  to 
have  a  very  small  gain  unbalance.  The  modulator  used  in  the  example  in  Section  F,  the 
ADL5375  from  Analog  Devices  [37],  specifies  an  I  and  Q  gain  unbalance  is  only  0.08 
dB. 


Figure  26.  Linear  Array  Pattern  Degraded  by  I  and  Q  Gain  Unbalance  (Solid). 

(a)  L  =  0.4  dB  and  (b)  L  =  1 .0  dB 
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Figure  27.  Relative  Average  Sidelobe  Level  {SIX')  vs.  I  and  Q  Gain 
Unbalance  L  Plotted  for  Four  Scanning  Angles 


Figure  28.  Relative  Peak  Sidelobe  Level  vs.  I  and  Q  Gain  Unbalance  L 

Plotted  for  Four  Scanning  Angles 
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When  compared  to  phase  noise  simulations  in  the  previous  section,  the  loss  due  to 
I  and  Q  gain  unbalance  is  small.  The  plot  in  Figure  29  shows  an  average  loss  of  0.002  dB 
when  the  gain  unbalance  is  greater  than  0.5  dB. 


Figure  29.  Array  Gain  Loss  Due  to  I  and  Q  Gain  Unbalance  L  Plotted  for 

Four  Scanning  Angles 


Beam  pointing  accuracy  is  more  affected  by  this  type  of  error  than  the  random 
phase  deviation.  The  average  beam  pointing  error  is  shown  in  Figure  30,  and  the  beam 
pointing  standard  deviation  relative  to  the  HPBW  is  shown  in  Figure  31.  From  the 

(BPE)  plot  one  can  see  that  the  pointing  error  is  more  significant  when  the  gain 

unbalance  is  higher.  The  pointing  error  standard  deviation  reaches  12%  with  0.5  dB  of 
gain  unbalance.  The  HPBW  as  a  function  of  gain  unbalance  is  shown  in  Figure  32. 
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Figure  30.  Relative  Average  Beam  Pointing  Error  { BPE }  vs.  I  and  Q  Gain 
Unbalance  L  Plotted  for  Four  Scanning  Angles 


Figure  3 1 .  Relative  Standard  Deviation  a vs.  I  and  Q  Gain  Unbalance  L 

Plotted  for  Four  Scanning  Angles 
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Figure  32.  Half-Power  Beamwidth  vs.  I  and  Q  Gain  Unbalance  L  Plotted  for 

Four  Scanning  Angles 


4.  Linear  Frequency  Drift 

a.  Simulation  Scenario 


For  linear  frequency  drift,  it  is  assumed  that  all  elements  have  an  independent  LO. 
If  no  other  phase  errors  are  considered,  the  instantaneous  frequency  (4.1)  becomes 


f(t)  =  fs+fd(t)  +  fa(t) 

=  /„+/rf(0 


(5.9) 


With  foffset=°  and/a( 0=0- 


For  this  simulation,  the  elements  are  assumed  to  have  independent  LOs  with 
instantaneous  frequency  normally  distributed  in  a  random  process,  with  PDF 

.23 


PDFfreq  = 


1 


exp 


( fzJn y 

2  °d 


(5.10) 


where 
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/  =  instantaneous  frequency / (/), 

f  =  nominal  frequency,  and 

<rd  =  frequency  drift  standard  deviation 

=  2jrDt2l 2.515. 

Recall  that =  2 nDt2. 

Using  this  definition  for  cr ^  leads  to  a  fj{t)  contribution  within  the  desired 
oscillator  accuracy  with  99%  confidence.  Furthermore,  when  t  =  0 ,  all  the  elements  are 
phase  synchronized  to  the  nominal  frequency /(/)  =  fn.  The  expression  for  the  array 
factor  accounting  only  for  frequency  drift  is  given  by 


AF(d,(f))  -  X  Wnejk'ron 


(5.11) 


where 


k*ron  =  k[uxn  +  vyn  +  wzn  ] 


(2  KUn  +  fd]/ c)[uxn  +  vyn  +  wzn]. 


Expanding  the  expression  (5.1 1),  we  get 


n= 1 


(5.12) 


where  kn  —  27rfnjc.  For  a  linear  array  along  the  x-axis  (5.12)  reduces  to 
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N 


AF(0,4>)  =  X  Wnejkn'r°n  exp 


f  j  lKfd(t^UXr  ^ 


n— 1 


V 


(5.13) 


J 


where  the  expression  in  parenthesis  is  the  phase  error  due  to  frequency  drift.  Notice  the 
expressions  (5.12)  and  (5.13)  are  actually  functions  of  time  because  the  linear  drift  is  time 
dependent.  An  important  observation  is  that  drift  contribution  to  phase  error  is  not  a 
random  process.  In  other  words,  independent  of  what  type  of  oscillator  is  being  used,  if 
they  are  set  in  free-run  or  hold-over  mode,  the  phase  error  magnitude  linearly  increases. 


b.  Simulation  Results 

Different  from  the  previous  scenarios,  for  frequency  drift  the  pattern  degradation 
is  illustrated  as  a  function  of  time.  The  importance  of  a  stable  local  oscillator  distribution 
system  to  continuously  synchronize  the  system  is  respectively  illustrated  in  this  section. 

In  Figure  33  and  Figure  34  the  gain  loss  and  average  sidelobe  level,  respectively, 
as  a  function  of  time  without  synchronization  of  elements  whose  LO  are  in  a  free-run 
mode  are  shown.  Notice  the  best  clock,  with  0.1  ppm  accuracy,  tolerates  a  longer  non¬ 
synchronization  interval,  around  0.1  s.  For  example,  it  is  necessary  to  employ  a  LO  with 
10-13  accuracy  if  an  array  requires  a  synchronized  period  of  24  hours  without  external 
reference.  The  cesium  frequency  standard  Microsemi  5071A  has  10-12  accuracy  [38]; 
moreover,  the  synchronization  requirement  is  more  demanding  and  critical  for  higher 
frequencies. 

In  Figure  35  an  array  pattern  simulation  assuming  an  array  employing  local 
oscillators  with  accuracy  of  0.1  ppm  or  better,  in  free-run  mode,  and  a  nominal  frequency 
of  300  MHz  is  shown.  An  array  using  a  LO  with  an  accuracy  of  100  ppm  was  simulated, 
and  the  results  are  illustrated  in  Figure  36.  This  scenario  is  more  realistic  because  100 
ppm  is  approximately  the  accuracy  of  available  COTS  VCO’s.  In  this  case  the  pattern  is 
severely  degraded  after  2.0  ms,  meaning  that  the  synchronization  interval  should  not  be 
longer  than  a  fraction  of  millisecond. 
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The  conclusion  extracted  from  the  Figure  35  and  Figure  36  pattern  plots  is  the 
necessity  of  periodic  phase  synchronization  in  order  to  avoid  pattern  degradation  due  to 
frequency  drift  errors. 


Figure  33.  Gain  Loss  Due  to  LO  in  Free-Run  as  a  Function  of  Time.  Different 
Colors  Represent  Different  LO  Accuracies,  Operating  at  300  MHz. 


s  30 


i  . . 1  i  i— r-r-m-r 

i  i  i  1 1 1 1 1| 

A 

*  p 

W * 

. 

■a— 0.1  ppm 
•••<#— 10  ppm 
a  -  1000  ppm 

_ jgj 

4 

4) 

4 

£ 

a 

1 

* 

A 

£ 

* 

i 

i 

_ ,  ,  ,  , 

0-0 

r-***;m 

M 

I 

10  10 
tkne  (s) 


10  10 


Figure  34.  Average  Sidelobe  Level  Due  to  LO  in  Free-Run  as  a  Function  of 
Time.  Different  Colors  Represent  Different  LO  Accuracies,  Operating 

at  300  MHz. 
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Figure  35.  Pattern  Degradation  Due  to  LO  in  Free-Run.  LO  Accuracies  of  0.1 
ppm  or  Better  (a)  t  =  1.8  fis  (b)  t=  1.9  ms  and  (c)  t  =  920  ms 


Figure  36.  Pattern  Degradation  Due  to  LO  in  Free-Run.  LO  Accuracies  of  100 
ppm  or  Better  (a)  t  =  1.8  fis  (b)  t  =  1.9  ms  and  (c)  t  =  920  ms 
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5. 


Element  Mislocation  Error 


a.  Simulation  Scenario 

In  Table  2,  the  simulation  settings  for  a  uniform  linear  array  of  50  identical 
elements  employed  in  the  element  dislocation  error  scenario,  which  has  a  different 
number  and  spacing  of  elements  compared  to  previous  simulations,  thereby  allowing  a 
larger  variation  in  element  mislocation,  is  summarized. 


Table  2.  Array  Characteristics  used  in  MATLAB  Simulations 


Number  of  elements 

50 

Elements  spacing 

0.4  wavelength 

Element  factor 

cos(0) 

Nominal  frequency 

300  MHz 

Tapering 

Taylor  (-35  dB  SLL,  Nbar  =  5) 

In  this  scenario,  a  single  LO  is  distributed  to  the  array  elements  in  a  series  fed 
manner  from  one  edge  of  the  array.  The  elements  are  randomly  mislocated  within  a  range 
so  that  two  adjacent  elements  do  not  overlap.  An  example  is  illustrated  in  Figure  37.  With 
a  series  feed  the  LO  distribution  errors  accumulate. 


Figure  37.  (a)  Ideal  Linear  Array,  (b)  Linear  Array  with  Elements  non- 

Uniformly  Distributed 
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In  an  ideal  scenario,  when  using  a  common  LO,  the  elements  can  compensate  for 
the  delay  between  element  locations  by  applying  a  correction  phase  for  the  scanning 
direction;  hence,  for  this  scenario  the  complex  weight  is  set  to  phase  compensate  based 
on  the  elements’  theoretical  positions.  Suppose  the  element  position  is  only  perturbed  in 
the  x  direction  so  yn  =  zn  =  0 ;  hence,  (3.20)  can  be  written  as 

rn  ~  ron  ~  Axn)x.  (5.14) 


The  correction  phase  in  the  scanning  direction  is 


Wn 


2jtxonsm(0s) 

2 


(5.15) 


The  distance  between  the  ideal  position  and  the  real  position  represents  a  phase  error  that 
is  expressed  as 

k  •  Ar«  =  2nPe  sin(0)  (5.16) 

where  P  is  the  position  error  and  is  defined  as  P  =  Ax^/2,  which  is  expressed  in 
fractions  of  wavelengths.  The  variable  P  is  a  uniform  random  variable  with  probability 
density  function 


PDF  = 


l 


MAX (Pf  )~  MIN(  P( ) 
0 


MIN{Pe )  <Pe<  MAX(Pe  ) 
otherwise 


(5.17) 


The  mislocation  of  elements  changes  the  array  factor,  which  is  expressed  as 


AF(6,<p )  =  X  Wnejk'f°»eJk'Af:"  ^ej2nsm^)P^  (5.18) 

n—\  n=l 

where  the  expression  in  parentheses  accounts  for  the  phase  error  due  to  element 
mislocation. 
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b. 


Simulation  Results 


In  Figure  38  we  see  the  pattern  disrupted  by  the  random  mislocation  of  the 
elements  and  additional  phase  noise.  On  the  top  plot  the  phase  noise  standard  deviation  is 
0.006  wavelengths,  whereas  the  standard  deviation  for  the  bottom  plot  is  0.058 
wavelengths. 

From  Figure  39  and  Figure  40,  we  see  there  is  significant  increase  in  average  and 
peak  sidelobes.  For  RMS  position  error  P(  >  0.03 A,  the  (SLL)  is  already  larger  than  the 
maximum  sidelobe  level  limited  by  Taylor’s  taper.  When  P  >  0.04 A,  the  peak  SLL  is  at 
least  10  dB  higher  than  the  no  error  scenario  and  is  worse  for  oblique  scanning  angles.  In 
Figure  41  the  array  gain  loss  due  to  element  mislocation  is  shown.  As  expected  the 
increase  in  SLL  increases  the  gain  loss;  however,  the  array  was  set  to  make  the  phase 
correction  in  scanning  direction,  which  reduces  the  loss  effect. 


Pattern  Angle,  6  (Degrees) 
(a) 


Figure  38.  Linear  Array  Aattern  Degraded  by  Element  Mislocation  (a)  RMS 
Mislocation  Error  0.006  Wavelength  and  (b)  RMS  Mislocation  Error 

0.058  Wavelength 
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Figure  39.  Relative  Average  Sidelobe  ( SLL)vs .  RMS  Element  Mislocation 
(Pe)  Expressed  in  Fraction  of  Wavelength  Plotted  for  Four  Different 

Scanning  Angles 


Figure  40.  Relative  Peak  Sidelobe  vs.  RMS  Element  Mislocation  (P ) 

Expressed  in  Fraction  of  Wavelength  Plotted  for  Four  Different 

Scanning  Angles 
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Figure  41.  Array  Gain  Loss  Due  and  RMS  Element  Mislocation  (P£)  Plotted 

for  Four  Scanning  Angles 


In  Figure  42  and  Figure  43  the  beam  pointing  error  average  and  standard 
deviation,  respectively,  are  shown.  The  latter  is  relative  to  the  half-power  beamwidth. 
The  beam  pointing  error  in  both  plots  increases  significantly  with  position  error.  For 
example,  the  CTgp^-  increased  seven  times  when  the  P£  changed  from  0.1  to  0.2  A  for 

0S  =  45° .  The  half-power  beamwidth  is  shown  in  Figure  44  and  does  not  show  any 
significant  change  due  to  positioning  error. 
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Figure  42.  Magnitude  of  Average  Beam  Pointing  Error  (BPEjv s.  RMS 
Element  Mislocation  (Pe)  Plotted  for  Four  Different  Scanning  Angles 


Figure  43.  Relative  Standard  Deviation  Beam  Pointing  Error  <JBpE  vs.  RMS 
Element  Mislocation  (Pe)  Plotted  for  Four  Different  Scanning  Angles 
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Figure  44.  Half -Power  Beamwidth  vs.  RMS  Element  Mislocation  (Pe) 
Plotted  for  Four  Different  Scanning  Angles 


E.  ARRAY  OF  SUBARRAYS  HARDWARE  ERROR  SIMULATIONS 

1.  Simulation  Scenario  for  Phase  Noise,  Gain  Unbalance  and  Frequency 
Drift  Error 

In  this  scenario  the  pattern  degradation  due  to  phase  error  for  an  array  of  linear 
subarrays  with  ten  identical  elements  per  subarray  was  simulated.  The  subarrays  are 
randomly  mislocated  within  four  and  five  wavelengths;  hence,  there  is  no  subarray 
overlap.  In  Table  3  and  Table  4,  the  array  characteristics  are  summarized. 
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Table  3.  Array  of  Subarrays  Characteristics  used  in  MATLAB  Simulations 


Number  of  subarrays 

20 

Subarray  spacing 

4-5  wavelengths  (uniform  random 
distributed) 

Tapering 

None 

Table  4.  Subarray  Characteristics  used  in  MATLAB  Simulations 


Number  of  elements 

10 

Elements  spacing 

0.4  wavelength 

Element  factor 

cos(0) 

Nominal  frequency 

300  MHz 

Subarray  tapering 

Taylor  (-20  dB  SLL,  Nbar  =  2) 

2.  Phase  Noise  Simulation  Results 

The  same  simulation  setting  used  for  the  linear  array  was  repeated  for  an  array  of 
subarrays;  i.e.,  each  element  has  independent  LOs  which  are  synchronized.  The  subarray 
pattern  with  added  phase  noise  deviation  q>  is  given  in  (3.8),  whereas  the  overall  array 
pattern  was  computed  using  (3.6). 

The  simulation  results  are  obtained  from  the  average  of  100  trials  for  four 
different  scanning  angles:  0°,  15°,  30°  and  45°. 

In  Figure  45(a)  the  pattern  distorted  by  phase  noise  with  RMS  of  0.05  rad  is 
shown,  and  in  the  bottom  plot,  Figure  45(b),  the  phase  noise  RMS  value  is  0.5  rad.  When 
compared  to  the  dotted  line,  which  is  the  pattern  without  noise  errors,  the  plot  in  Figure 
45(a)  is  almost  indistinguishable,  whereas  for  the  plot  in  Figure  45(b),  the  increase  in 
sidelobe  level  is  noticeable.  This  is  confirmed  by  the  relative  average  and  peak  sidelobes 
plots  in  Figure  46  and  Figure  47.  Note  the  sidelobes  for  array  of  subarrays  are  usually 
higher  than  a  linear  array  because  of  the  grating  lobes  caused  by  the  subarray  spacing. 
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In  Figure  48  the  gain  loss  due  to  random  phase  deviation  is  illustrated.  For  phase 
noise  with  RMS  of  one  radian  (a  very  large  value),  the  gain  loss  is  approximately  3  dB. 
Overall,  the  values  for  gain  loss  are  similar  to  those  obtained  in  the  linear  array 
simulation. 


(b) 


Figure  45  Linear  Array  of  Subarrays  Pattern  Degraded  by  Phase  Noise 
(blue),  the  Dotted  Line  is  the  Array  Without  Phase  Noise  Error 
(a)  (pRK1K  =  0.05  rad  and  (b)  (pRMX  =  0.5  rad 


Figure  46.  Relative  Average  Sidelobe  (SLL  ')vs.  RMS  Phase  Deviation  (pj^s 
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Plotted  for  Four  Different  Scanning  Angles 


Figure  47.  Relative  Peak  Sidelobe  Level  vs.  RMS  Phase  Deviation 
Plotted  for  Four  Different  Scanning  Angles 


Figure  48.  Array  Gain  vs.  RMS  Phase  Deviation  Plotted  for  Four 

Scanning  Angles 


In  Figure  49  the  average  beam  pointing  error  is  illustrated  and  shows  a  small 
variation,  on  the  order  of  couple  of  millidegrees,  caused  by  phase  noise.  Additionally,  in 
Figure  50  the  relative  beam  pointing  error  standard  deviation,  which  is  small  compared  to 
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the  result  illustrated  in  Figure  23  for  a  uniform  array,  is  shown.  In  the  current  simulation, 
the  half-power  beamwidth  is  smaller  because  the  scenario  used  a  larger  array  and  more 
elements. 


Figure  49.  Average  Beam  Pointing  Error  I^BPE^v s.  RMS  Phase  Deviation 
<PrMs  for  Four  Different  Scanning  Angles 


Figure  50.  Relative  Standard  Deviation  Beam  Pointing  Error  <JppE  vs.  RMS 
Phase  Deviation  Plotted  for  Four  Different  Scanning  Angles 
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Figure  51.  Half-Power  Beam  width  vs.  RMS  Phase  Deviation  <Prms  Plotted 

for  Four  Different  Scanning  Angles 

3.  I  and  Q  Gain  Unbalance  Simulation  Results 

In  this  simulation  the  pattern  degradation  caused  by  I  and  Q  gain  unbalance  was 
obtained;  hence  the  subarray  pattern  accounting  for  these  errors  is  given  in  (5.8),  and  the 
overall  array  pattern  was  computed  using  (3.6). 

Similar  to  the  previous  case,  the  results  in  this  section  are  the  average  of  100  trials 
for  each  scanning  angle,  i.e.  0°,  15°,  30°  and  45°. 

From  the  array  patterns  illustrated  in  Figure  52  one  can  conclude  that  this  array 
configuration  is  not  affected  by  small  gain  errors  since  Figure  52(a)  shows  a  pattern  for 
an  array  with  gain  errors  up  to  1  dB.  Although  the  pattern  shown  in  Figure  52(b)  presents 
high  sidelobes,  to  get  this  result  the  simulation  considered  gain  unbalance  up  to  5  dB, 
which  is  unlikely  to  happen  except  for  the  major  hardware  failures.  For  instance,  the 
specification  of  modulator  ADL5375  from  Analog  Devices  [37]  for  I  and  Q  gain 
unbalance  is  only  0.08  dB. 
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Figure  52.  Linear  Array  of  Subarrays  Pattern  Degraded  Because  of 
Unbalanced  Gain  (blue),  and  the  Dotted  Line  is  the  Array  without 
Errors  (a)  L=  1.0  dB  and  (b)  L  =  5.0  dB 


In  Figure  53  and  Figure  54  the  relative  average  and  peak  sidelobes  plots, 
respectively,  are  shown.  The  sidelobe  level  increase  is  significant  only  for  large  gain 
unbalance.  This  confirms  the  previous  statement  that  this  type  of  error  does  not 
significantly  affect  the  sidelobe  level. 

Regarding  the  gain  loss  results,  when  we  compare  the  plots  in  Figure  29  and 
Figure  55,  the  loss  for  a  linear  array  is  slightly  better  than  that  for  an  array  of  subarrays. 
On  the  other  hand,  this  simulation  showed  that  the  beam  pointing  performance,  both 
average  and  relative  standard  deviation,  for  an  array  of  subarrays  is  better  than  the  linear 
array  when  the  plots  from  Figure  30  and  Figure  31  are  compared  with  those  in  Figure  56 
and  Figure  57,  respectively.  The  better  beam  pointing  performance  is  a  consequence  of  a 
larger  array  when  the  I  and  Q  gain  errors  are  normally  distributed.  If  this  assumption  is 
true,  and  these  errors  are  spatially  uncorrelated,  then  the  average  approaches  zero  at  the 
scanning  angle.  Finally,  the  graph  in  Figure  58  shows  the  half-power  beamwidth  as  a 
function  of  gain  unbalance,  which  obviously  does  not  affect  this  characteristic. 
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Peak  SLL  MAX(SLL)  (dB) 


-23 


Figure  53.  Relative  Average  Sidelobe  (iS'ZX)lncrease  Due  to  I  and  Q  Gain 
Unbalance  Plotted  for  Four  Different  Scanning  Angles 


Figure  54.  Relative  Peak  Sidelobe  Level  Increase  Due  to  I  and  Q  Gain 
Unbalance  Plotted  for  Four  Different  Scanning  Angles 
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Figure  55.  Array  Gain  Loss  Due  due  to  I  and  Q  Gain  Unbalance  Plotted  for 

Four  Scanning  Angles 


Figure  56.  Average  Beam  Pointing  Error  (BPE)  vs.  I  and  Q  Gain  Unbalance 

Plotted  for  Four  Different  Scanning  Angles 
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Figure  57.  Relative  Standard  Deviation  Beam  Pointing  Error  <rgpE  vs.  I  and 
Q  Gain  Unbalance  Plotted  for  Four  Different  Scanning  Angles 


/  and  Q  gain  unbalance  (dB) 


Figure  58.  Half-Power  Beamwidth  vs.  I  And  Q  Gain  Unbalance  Plotted  for 

Four  Different  Scanning  Angles 
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F. 


TYPICAL  SCENARIO 


The  objective  of  this  section  is  to  show  an  example  of  the  impact  of  phase  noise 
and  gain  unbalance  on  the  array  pattern  using  the  information  obtained  from  the  data 
sheet  of  the  COTS  components. 

Suppose  a  surface  ship  radar  is  operating  at  S-band,  i.e.,  3500  MHz,  employing  a 
uniform  digital  phased  array  design  with  one  hundred  elements,  which  are  individually 
affected  by  phase  errors.  For  this  analysis  the  system  bandwidth  is  10  MHz.  Shown  in 
Figure  59  is  a  block  diagram  of  a  T/R  module  with  non-ideal  components,  which  are 
listed  in  Table  5. 


Figure  59.  T/R  Module  Block  Diagram  Used  as  Example  of  Digital  Phased 

Array  with  Non-Ideal  Components 
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Table  5.  Example  T/R  Module  Components  and  Function 


Function 

Manuf. 

Model 

Reference 

Oven  controlled 
crystal  oscillator 

Synchronization  reference 
sent  by  DBF  to  elements 

Wenzel 

501-26636 

[39] 

Phased  locked  loop 
(VCO) 

LO  reference  to  digital 
components,  modulator 
and  demodulator. 

Hittite 

HMC834 

LP6GE 

[40] 

Digital-to-Analog 
Converter  and 
Waveform  Generator 

Digital  waveform 
generator. 

Analog 

Devices 

AD9102 

[41] 

Analog-to-Digital 

Converter 

Digitalization  of  received 
baseband  signal. 

Analog 

Devices 

AD9642 

[42] 

Quadrature 

modulator 

Modulate  and  upconvert 
the  waveform  using  LO  as 
carrier  multiplier. 

Analog 

Devices 

ADL5375 

[37] 

Quadrature 

demodulator 

Demodulate  the  received 

RF  signal  using  LO  as 
carrier  multiplier. 

Analog 

Devices 

ADL5380 

[43] 

1.  Calculating  Phase  Noise 

In  order  to  convert  the  phase  noise  specified  in  oscillators’  datasheets  to  random 
phase  (p,  also  defined  as  jitter  in  [44],  the  phase  noise  PSD  has  to  be  integrated  over  the 
signal  bandwidth,  which  is  equivalent  to  computing  the  area  under  the  curve  shown  in 
Figure  60: 


(Phase  Noise  Power)  =  Arca^  +  Area^  +  Areac  +  Area^  +  Areaf.  (5.19) 

The  trapezoidal  area  marked  in  Figure  60  is  an  area  on  a  log-log  plot,  which  has  x  and  y 
axis  in  logarithmic  scale.  The  total  random  phase  deviation  is  computed  using  the  square- 
root  sum  of  the  individual  phase  noise  contributors 


<Prms  ~ 


£-lfinV10  i  N 

21  10  * 

i=  1  V 


m,  + 1 

Ji 


^/w2+l  /7W2+ 1 

Jill  ~Ji 


(5.20) 
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where  the  frequencies  f  are  the  corner  frequencies  between  phase  noise  regions,  K  is 
the  number  of  phase  noise  regions  and  i  are  the  curve  slopes  defined  as 


y/,+i)-y  fd 

m‘  logio(/i+i)-logio(/i) 


(5.21) 


Figure  60.  Phase  Noise  Frequency-Domain  Representation  of  A  MCXO 
Model  Vectron  MX-503.  The  Colored  Area  under  the  Curve 
Represents  the  Phase  Noise  Power 

Adapted  from  [1]:  W.  Kester,  (2008,  Dec.).  “Converting  Oscillator  Phase  Noise  to  Time 
Jitter,”  2008.  [Online],  Available:  http://www.analog.com/media/en/training-seminars/ 
tutorials/MT -008 .pdf. 


The  upper  limit  for  integration  (number  of  phase  noise  regions)  is  determined  by 
signal  single-sided  bandwidth,  whereas  the  lower  limit  is  found  on  the  LO  data  sheet.  In 
practice,  data  sheets  usually  provide  phase  noise  values  from  100  Hz  to  1  MHz  of 
frequency  offset,  neglecting  the  phase  noise  close  to  the  carrier,  because  the  majority  of 
the  jitter  is  due  to  the  WPM  noise.  Alternatively,  there  are  online  calculators  and 
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spreadsheets  that  convert  the  phase  noise  to  RMS  (Pj^s  ra<hans  or  jitter  in  seconds, 
like  the  one  found  in  [45]. 

In  order  to  calculate  the  cumulative  effect  of  phase  noise  for  two  or  more 
components,  the  noise  power  of  the  individual  component  contributions  is  summed 


where 


N  K  1 

<PRMS  =  P-'Yj  2 

J  I  *  I 


s,  /to  ^ 

10  1 


f- 

Jy 


mi  mu  + 1 


y 


fmij 1 1  _  ,mij+l 

Jij+l  Jij 


(5.22) 


fy  -  the  comer  frequencies  between  phase  noise  regions 
of  the  7 th component  in  the  cascade, 

nty  =  i“  phase  noise  region  curve  slope  of  the  jtb  component 

in  the  cascade,  and 

N  =  number  of  components  in  cascade. 

2.  Crystal  Oscillator  Wenzel  501-26636 

In  the  proposed  scenario  the  XO  Wenzel  501-26636  is  the  reference  transmitted 
from  the  DBF  to  the  T/R  modules  for  synchronization.  From  the  manufacturer  data  sheet 
the  XO  phase  noise  is  listed  in  Table  6. 

Using  (5.20),  we  get  the  RMS  value  for  random  phase  deviation  at  320  MHz: 


^,/?MS'|y-320MHz  1-135x10  rad. 


(5.23) 
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Table  6.  XO  Wenzel  501-26636  Phase  Noise  Specification 


100  Hz 

-118  dBc/Hz 

1  kHz 

-141  dBc/Hz 

10  kHz 

-151  dBc/Hz 

100  kHz 

-152  dBc/Hz 

3.  Voltage  Controlled  Oscillator  Hittite  HMC834LP6GE 

The  reference  signal  from  the  XO  Wenzel  501-26636  is  used  for  frequency 
synthesis  using  VCO  Hittite  HMC834LP6GE,  whose  phase  noise  specification  is  shown 
in  Table  7  and  Table  8. 


Table  7.  VCO  Hittite  HMC834LP6GE  Phase  Noise  Specification  at  4  GHz 


10  kHz 

-78  dBc/Hz 

100  kHz 

-108  dBc/Hz 

1  MHz 

-134  dBc/Hz 

10  MHz 

-156  dBc/Hz 

100  MHz 

-171  dBc/Hz 

At  4  GHz,  the  RMS  value  for  random  phase  deviation  is 

<PRMs\ /=  4GHz  =  0  0126  (5  24) 
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Table  8.  VCO  Hittite  HMC834LP6GE  Phase  Noise  Specification  at  100  MHz 


10  kHz 

-111  dBc/Hz 

100  kHz 

-141  dBc/Hz 

1  MHz 

-163.5  dBc/Hz 

10  MHz 

-170  dBc/Hz 

100  MHz 

-173  dBc/Hz 

At  100  MHz,  the  RMS  value  for  random  phase  deviation  is 

<Pkms\ y-iooMHz  =  2‘8312  x  10^  rad  (5-25> 

Adding  the  contribution  from  reference  XO,  we  get  the  LOs  reference  signal  for  the 
modulator  and  demodulator  with  a  random  phase  deviation  of 

PRMSlo  =  yj(<PRMs\f  320MHz)2 +(^«7  4GHz  ^ 

=7(1-135  x  10^)2  +(0.0126)2  (5.26) 

=  0.0126  rad. 

The  clock  reference  signal  for  the  DDS  and  A/D  converter  has  a  random  phase  deviation 
of 

<PRMS.it  =  ^PKMs\ /=320MHz)2  7=100^ 

=7(1-135  x  10^)2  +  (2.8312  x  10^)2 

=  3.0501x10"*  rad. 
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(5.27) 


Notice  for  low  frequencies  the  phase  noise  contribution  is  negligible;  on  the  other 
hand,  for  high  frequencies,  the  phase  noise  increases  significantly.  In  Table  9,  the 
expected  pattern  alteration  due  to  phase  noise  is  summarized. 


Table  9.  Pattern  Degradation  Due  to  LO  Phase  Noise 


os  =  v° 

e,  =  i5° 

*,  =  30° 

*,=45° 

Average  SLL  increase  (dB) 

0.4 

0.42 

0.39 

0.3 

Peak  SLL  increase  (dB) 

0.72 

0.64 

0.7 

0.71 

Loss  (dB) 

0.05 

0.05 

0.05 

0.05 

Relative  beam  pointing  error 
standard  deviation 

0.0035 

0.0035 

0.0039 

0.0044 

4.  Phase  Errors  Due  to  I  and  Q  Gain  Unbalance 

For  the  proposed  scenario,  the  modulator  ADL5375  and  the  demodulator 
ADL5380  have  the  I  and  Q  gain  unbalance  listed  in  Table  10.  Additionally,  the  DDS 
(AD9102)  and  A/D  converter  (AD9642)  contribute  with  a  gain  error  of  1%  and  2%, 
respectively.  This  is  equivalent  to  an  additional  gain  error  of  0.04  dB  in  transmission  and 
0.08  dB  on  reception. 


Table  10.  Modulator,  Demodulator  Gain  Unbalance  at  3500  MHz,  and  DDS  and 

A/D  Converter  Gain  Error 


Modulator  ADL5375  (dB) 

0.08 

Demodulator  ADL5380  (dB) 

0.14 

DDS  AD9102  (dB) 

0.04 

A/D  converter  AD9642  (dB) 

0.08 
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In  Table  11  and  Table  12,  the  pattern  degradation  due  to  transmission  and 
reception  gain  errors,  respectively,  are  summarized. 


Table  1 1 .  Pattern  Degradation  Due  to  Transmission  I  and  Q  Gain  Unbalance 


es  =  o° 

o 

IT) 

II 

CD 

=  30° 

0^=45° 

Average  SLL  increase  (dB) 

0.03 

0.03 

0.02 

0.03 

Peak  SLL  increase  (dB) 

0.06 

0.1 

0.14 

0.11 

Loss  (dB) 

3x1 0'5 

6x1 0'5 

3xl0'5 

2x1  O'4 

Relative  beam  pointing  error 
standard  deviation 

0.0316 

0.0316 

0.0316 

0.0316 

Table  12.  Pattern  Degradation  Due  to  Transmission  I  and  Q  Gain  Unbalance 


o 

O 

II 

<£> 

e,  =  i5° 

=  30° 

0S  =  45° 

Average  SLL  increase 

0.11  (dB) 

0.12  (dB) 

0.1  (dB) 

0.07  (dB) 

Peak  SLL  increase 

0.18  (dB) 

0.32  (dB) 

0.3  (dB) 

0.3  (dB) 

Loss 

5x1 0'4  (dB) 

lxlO"4  (dB) 

2x1 0~4  (dB) 

0  (dB) 

Relative  beam  pointing  error 
standard  deviation 

0.0608 

0.0608 

0.0608 

0.0608 

5.  Overall  Pattern  Degradation 

In  Table  13,  the  overall  pattern  degradation  due  to  phase  noise  and  gain  errors, 
accounting  for  transmission  and  reception,  is  summarized.  Although  the  components  used 
in  this  example  have  small  errors  individually,  the  pattern  degradation  due  to  the 
combined  effect  of  phase  noise  and  gain  error  is  not  negligible. 
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Table  13.  Overall  Pattern  Degradation 


o 

o 

II 

CD 

0S  =  15° 

05=3O° 

05  =  45° 

Average  SLL  increase  (dB) 

0.94 

0.99 

0.9 

0.7 

Peak  SLL  increase  (dB) 

1.68 

1.7 

1.84 

1.83 

Loss  (dB) 

0.1 

0.1 

0.1 

0.1 

Relative  beam  pointing  error 
standard  deviation 

0.0687 

0.0687 

0.0687 

0.0687 

G.  SUMMARY 

In  this  chapter,  the  array  pattern  degradation  was  analyzed  using  results  generated 
from  equations  programed  in  MATLAB  to  simulate  and  compute  the  pattern 
characteristics:  average  sidelobe  level,  peak  sidelobe  level,  gain  loss,  average  and  relative 
standard  deviation  for  beam  pointing  error,  and  HPBW. 

Phase  noise  and  I  and  Q  gain  unbalance  were  simulated  for  two  different 
scenarios:  a  uniform  linear  array  and  a  linear  array  of  subarrays.  In  both  cases,  the  results 
show  a  significant  increase  in  sidelobe  levels  but  small  gain  loss.  For  the  uniform  linear 
array  scenario,  the  beam  pointing  error  is  larger  than  for  a  linear  array  of  subarrays 
because  of  the  difference  in  size  between  the  arrays. 

For  phase  errors  caused  by  element  mislocation,  only  the  uniform  linear  array 
scenario  was  simulated.  The  results  show  that  small  deviations  can  produce  significant 
degradation  in  sidelobe  level  and  beam  pointing  error;  moreover,  the  pattern 
characteristics  are  worse  for  oblique  angles.  The  gain  loss  is  not  significantly  affected. 

In  Section  D.4,  the  relationship  of  the  non-synchronization  interval  to  the  non¬ 
ideal  oscillator  accuracy  was  addressed.  For  this  simulation,  the  oscillator  was  modeled  to 
have  linear  frequency  drift.  From  the  simulation  results,  one  can  conclude  that  the  pattern 
resilience  depends  on  the  oscillator  accuracy  and  the  system  nominal  frequency. 
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In  Section  F,  a  typical  scenario  analysis  for  a  uniform  linear  array  constituted  of 
COTS  components  was  presented.  In  the  next  chapter,  a  summary,  the  conclusions  of  this 
research,  and  suggestions  for  future  work  are  presented. 
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VI.  CONCLUSION  AND  FUTURE  WORK 


The  primary  objective  of  this  research  was  to  analyze  phased  array  digital 
beamforming  degradation  under  the  influence  of  phase  errors  caused  by  phase  noise,  7 
and  Q  gain  unbalance  and  element  mislocation.  Parametric  simulations  were  programmed 
in  MATLAB  to  perform  Monte  Carlo  simulations  of  hardware  errors  using  equations  and 
models  studied  in  the  literature  review.  The  simulation  results  can  be  used  as  reference  to 
establish  acceptable  tolerance  limits  for  phase  errors. 

In  the  first  four  chapters  of  this  thesis  the  literature  review  of  topics  related  to 
digital  phased  array  antennas  was  presented.  The  literature  review  focused  on  the 
technologies  and  techniques  that  affect  the  coherency  among  the  array  elements,  which  is 
vital  for  beamforming.  Appropriate  error  models  were  developed  and  implemented  in 
MATLAB. 

Analyzing  the  T/R  module  architecture,  two  main  sources  for  phase  errors  were 
identified:  phase  noise  and  gain  unbalance  between  7  and  Q  branches.  Every  digital  and 
analog  component  that  requires  a  reference  signal  is  affected  by  local  oscillator  phase 
noise.  Moreover,  gain  unbalance  between  7  and  Q  branches  can  be  a  source  of  phase 
errors  in  quadrature  modulators  and  demodulators. 

In  order  to  build  the  simulation  scenarios,  modifications  to  the  standard  digital 
array  pattern  equations  were  proposed  to  account  for  phase  noise,  7  and  Q  gain  unbalance 
and  element  mislocation. 

Additionally,  the  mathematical  models  for  instantaneous  phase  and  phase  noise 
were  presented  and  were  used  to  model  the  local  oscillator  used  in  the  simulations.  The 
phase  noise  model  was  based  on  the  power-law  model,  which  is  an  approximation 
frequently  used  for  oscillators. 
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A.  CONCLUSION 

The  array  pattern  degradation  was  observed  from  the  analysis  of  the  sidelobe  level 
(average  and  peak),  gain  loss  and  beam  pointing  error  (average  and  standard  deviation). 
Phase  noise,  I  and  Q  gain  unbalance,  element  mislocation  and  frequency  linear  drift  were 
the  phase  error  sources  studied  in  these  experiments. 

Initially,  the  established  simulation  scenario  consisted  of  100  elements  uniformly 
distributed.  The  element  separation  distance  was  0.2  wavelengths.  For  phase  noise  and  I 
and  Q  gain  unbalance,  there  was  a  noticeable  increase  in  sidelobe  levels  and  beam 
pointing  error  proportional  to  the  phase  error.  On  the  other  hand,  the  gain  loss  was  not 
significant  in  both  cases. 

Additionally,  for  the  same  scenario,  the  effect  of  linear  frequency  drift  on  the 
array  pattern  was  simulated.  From  the  simulation  results,  we  see  a  relationship  between 
pattern  distortion,  oscillator  accuracy  and  system  nominal  frequency.  If  a  less  accurate 
device  replaces  an  oscillator  that  has  less  frequency  drift,  the  synchronization  interval 
must  decrease  proportionally  to  maintain  the  pattern  characteristics.  In  a  similar  fashion, 
higher  frequency  systems  must  be  synchronized  more  often. 

For  element  mislocation  simulations,  the  scenario  was  slightly  different  (50 
elements,  uniformly  distributed  with  separation  distance  between  elements  of  0.4 
wavelength);  however,  the  observed  effects  were  the  same.  That  is,  a  significant  increase 
of  sidelobe  level  and  beam  pointing  error  resulted.  Additionally,  it  was  observed  that  for 
more  oblique  scanning  angles,  the  degradation  of  these  pattern  characteristics  was  more 
significant. 

The  experiments  for  phase  noise  and  I  and  Q  gain  unbalance  were  repeated  for  an 
array  of  subarrays  scenario,  with  20  subarrays  randomly  spaced  from  one  another  within 
four  and  five  wavelengths.  Each  subarray  was  set  to  have  ten  elements.  In  both  cases,  the 
simulation  results  showed  a  significant  increase  in  sidelobe  levels  due  to  grating  lobes  in 
addition  to  the  effects  of  the  simulated  phase  errors  (phase  noise  and  gain  unbalance),  but 
in  contrast  a  small  value  of  array  gain  loss  was  found.  Compared  to  the  uniform  linear 
array  scenario,  this  experiment  obtained  better  performance  for  beam  pointing  error 
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parameters,  which  is  attributed  to  the  size  of  the  array  and  the  assumption  of  spatially 
uncorrelated  errors  among  the  elements. 

In  addition  to  the  error  sources  examined  in  this  study,  there  exist  other  factors 
related  to  phase  errors  that  were  not  considered  in  this  research: 

1.  Digital  components  of  the  array  introduce  quantization  errors 

2.  System  thermal  noise  was  neglected 

3.  Transmission  and  reception  coupling 

In  Chapter  V,  Section  F,  a  typical  radar  scenario  was  analyzed  using  information 
from  COTS  components’  specifications.  Among  the  array  pattern  characteristics,  the 
sidelobe  level  was  most  affected  by  phase  errors  (average  increase  of  1  dB),  but  the 
significance  of  the  degradation  depends  upon  the  system  requirements.  The  beam 
pointing  error  was  also  significant  in  this  example,  but  it  decreases  if  the  size  of  the  array 
increases.  Overall  gain  is  not  significantly  affected  by  the  phase  errors  considered  in  this 
research.  Array  designers  can  use  a  similar  analysis  to  establish  phase  error  limits  based 
on  the  desired  pattern  characteristics. 

B.  FUTURE  WORK 

1.  Phase  Noise  Analysis  for  EHF  Systems 

More  recently,  the  demand  for  bandwidth  pushed  the  research  for  radar  systems  to 
higher  frequencies,  such  as  EHF  or  millimetric  waves,  which  requires  high  performance 
components;  therefore,  these  systems  require  multiple  upconversion  and  downconversion 
stages,  which  have  an  impact  on  phase  noise. 

Moreover,  each  of  these  stages  require  non-linear  devices  (mixers)  to  perform 
frequency  multiplication,  increasing  the  probability  of  IM  products  being  generated. 
These  spurious  frequencies  increase  the  correlated  phase  noise,  which  can  degrade  the 
array  pattern  even  more. 
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2.  LO  Distribution  Techniques 

The  transmission  of  the  LO  signal  through  a  channel  and  its  recovery  using  non¬ 
ideal  devices  introduce  noise  that  is  transferred  to  other  components  in  form  of  phase 
noise.  Moreover,  the  synchronization  process  also  adds  phase  errors  due  to  discrete  phase 
shifters.  It  is  important  to  model  the  impact  of  these  errors  on  phase  noise  and  account  for 
them  in  the  array  pattern. 

Alternatively,  the  array  elements  can  be  locally  synchronized  by  a  LO  distribution 
system  divided  into  levels.  The  first  level  elements  are  synchronized  to  a  LO,  which  is 
linked  to  the  reference  of  LO  in  level  two.  These  second  level  LOs  are  directly  connected 
to  the  DBF.  Depending  on  the  size  of  the  array,  we  add  more  levels  of  LO;  consequently, 
the  random  errors  are  also  added  in  a  chain  effect. 


96 


APPENDIX.  LOCAL  OSCILLATORS 


Some  of  the  typical  oscillator  architectures  and  technology  and  examples  of 
available  COTS  devices  are  described  in  this  section. 


Crystal  Oscillators 


a. 


In  reference  [46],  the  author  describes  the  construction  and  operation  of  four  types 
of  crystal  oscillators. 

(1)  Crystal  Oscillators 

A  crystal  oscillator  (XO)  usually  has  quartz  as  a  resonance  element.  Quartz  is  a 
monocrystal  of  silicon  dioxide  (SiOa),  which  is  cut  and  packetized  to  supply  a  constant 
and  stable  resonance  frequency.  To  differentiate  from  other  XOs,  they  are  also  called 
packaged  crystal  oscillators  (PXOs)  [46]. 

The  limit  for  resonance  frequency  depends  on  the  crystal  mode  of  oscillation.  In 
Figure  61  (a)  and  (b)  are  illustrated  the  two  low  frequency  modes,  flexural  and  torsional, 
commonly  found  in  tuning  fork  resonators.  For  medium  frequencies  the  contour  mode  is 
usually  encountered,  exemplified  in  Figure  61  (c);  to  obtain  higher  frequencies  the 
thickness  mode  is  used,  which  is  illustrated  in  Figure  61  (d).  For  the  latter  mode, 
resonance  frequency  is  inversely  proportional  to  the  crystal  thickness,  so  the  cut  of  the 
crystal  also  influences  in  its  performance  [17]. 


(a)  flexural  mode 


tc)  contour  modes 


(b)  torsional  mode 


(d)  thickness  modes 


Figure  61.  The  Four  Possible  Modes  of  Oscillation  for  Quartz  Crystals 

Source  [17]:  E.  Vittoz,  Low-Power  Crystal  and  MEMS  Oscillators:  The  Experience  of 
Watch  Developments.  Lausanne:  Springer  Science+Business  Media  B.V.,  2010,  p.16. 
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(2)  Temperature  Compensated  Crystal  Oscillators 

Temperature  variations  might  change  the  resonator  characteristics,  degrading  the 
oscillator  accuracy.  Temperature  compensated  crystal  oscillators  (TCXO)  use  a 
temperature  sensor  signal  to  control  a  variable  capacitor  in  the  oscillator  network; 
therefore,  the  reactance  is  adjusted  despite  the  temperature  variation  [46],  According  to 
[46]  a  TCXO  is  20  times  more  stable  than  a  XO. 

(3)  Oven  Controlled  Crystal  Oscillators 

Oven  controlled  crystal  oscillators  (OCXO)  have  an  oven,  which  is  an  enclosed 
structure,  to  stabilize  the  sensitive  elements  of  the  oscillator,  such  as  the  crystal.  The  oven 
is  adjusted  to  a  temperature  where  the  crystal  is  stable.  The  OCXO  is  1000  times  more 
stable  than  XO  [46]. 

(4)  Microcomputer  Compensated  Crystal  Oscillators 

The  microcomputer  compensated  crystal  oscillators  (MCXO)  use  a 
microcomputer  to  measure  and  compensate  the  temperature  effects  over  the  resonator. 
There  are  two  types  of  MCXO,  one  using  pulse  deletion  and  other  using  frequency 
summation.  The  former  method  has  a  resonator  that  is  tuned  slightly  above  the  desired 
frequency;  then  pulses  are  eliminated  to  compensate  the  temperature  effect.  For 
frequency  summation,  the  microcomputer  uses  a  DDS  to  generate  a  frequency 
compensation  to  be  added  to  resonator  output  frequency;  the  result  is  a  constant 
frequency  [46]. 

b.  Atomic  Standards  and  Chip  Scale  Atomic  Clock 

All  atomic  frequency  standards  available  on  the  market  take  advantage  of 
hydrogen-like  atom  gases,  which  have  one  “unpaired  electron  in  a  symmetric  orbit”  [46], 
In  short,  the  physical  concept  behind  the  atomic  resonator  consists  of  a  beam  radiation  at 
the  frequency  that  causes  an  electron  excitation  on  a  hydrogen-like  gas,  which  forces  the 
atoms  to  change  their  energy  level.  Based  on  the  element  characteristics,  some  different 
techniques  are  used  to  extract  the  energy  state  frequency  information  and  are  illustrated 
in  the  next  sections. 
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(5)  Cesium  Beam  Frequency  Standard 

The  cesium  beam  frequency  standard  (CsFS)  is  illustrated  in  Figure  62.  The  CsFS 
has  a  cesium  gas  that  is  accelerated,  forming  an  atomic  beam,  which  is  submitted  to  a 
magnetic  field  that  selects  the  lower  energy  state  atoms.  These  atoms  form  a  beam  that  is 
directed  to  a  microwave  cavity.  Once  the  correct  frequency  is  transmitted  inside  the 
cavity,  the  atoms  are  stimulated  to  change  energy  state.  A  probe  is  used  to  transmit  the 
signal  in  the  cavity.  A  high-energy  state  cesium  beam  is  formed,  and  then,  it  is  submitted 
to  another  magnetic  field.  When  the  beam  hits  the  detector,  the  system  knows  the 
frequency  of  9,192,631,770  Hz  is  in  the  cavity  [47], 
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Figure  62.  CsFS  Operation.  Left  the  Atomic  State  Selection.  Right  Beam 

Detection 

Source  [46]:  J.  R.  Vig,  (2015,  NovlO)  Quartz  Crystal  Resonators  and  Oscillators  -  For 
Frequency  Control  and  Timing  Applications — A  Tutorial.  [Online].  Available: 
http://www.ieee-uffc.org/frequency-control/learning-vig-tut.asp. 


(6)  Rubidium  Frequency  Standard 

Rubidium  frequency  standard  (RbFS)  is  another  type  of  atomic  standard  that  uses 
characteristics  of  the  rubidium  (Rb)  atom  that  absorbs  a  specific  radiation  frequency  to 
change  state.  In  Figure  63  is  illustrated  a  RbFS  block  diagram.  The  87Rb  lamp  radiates 
light  that  can  be  absorbed  by  Rb-87  buffer  gas  with  high-energy  state.  When  the  radiation 
in  the  cavity  has  the  exact  frequency  of  6.834,682,608  GHz,  the  Rb  changes  from  low 
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energy  state  to  high.  Consequently,  the  light  is  partially  absorbed  by  the  gas  and  power 
on  the  photo  cell  drops.  A  quartz  VCXO  is  used  to  control  the  frequency  in  the  cavity.  To 
work  properly,  the  light  must  not  have  the  795  nm  wavelength  component  because  it 
prevents  the  Rb  atom  from  changing  state;  for  this  reason  a  filter  of  Rb  is  used  [46]. 
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Figure  63.  Rbfs  Operation  Block  Diagram 

Adapted  from  [46]  J.  R.  Vi g,  (2015,  NovlO)  Quartz  Crystal  Resonators  and  Oscillators  - 
For  Frequency  Control  and  Timing  Applications  -  A  Tutorial.  [Online].  Available: 
http://www.ieee-uffc.org/frequency-control/learning-vig-tut.asp. 


(7)  Hydrogen  Maser  Frequency  Standard 

In  a  Hydrogen  maser  standard  (H-maser),  a  beam  of  hydrogen  atoms  is  formed  by 
electrical  discharge  on  the  gas.  A  magnetic  field  is  used  to  select  the  beam  of  atoms  that 
has  a  certain  energy  state,  and  it  is  directed  to  a  microwave  cavity.  The  cavity  is  excited 
with  a  FM  signal  at  1,420,405,751.77  Hz,  which  induces  the  atoms  to  transition. 
Consequently,  the  cavity  transmits  an  AM  signal  through  a  probe,  and  the  desired 
frequency  is  extracted  from  the  signal  envelope  [48]. 

(8)  Chip  Scale  Atomic  Clock 

The  great  challenge  for  the  development  of  a  chip  scale  atomic  clock  (CSAC)  is 
the  fact  that  all  conventional  technologies  use  a  microwave  cavity,  limiting  the  dimension 
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reduction.  The  solution  to  this  is  the  coherent  population  trapping  (CPT)  technique.  In 
Figure  64  is  illustrated  the  block  diagram  of  a  CSAC. 


The  vertical-cavity  surface  emitting  lasers  (VCSELs)  transmit  two  laser  beams, 
whose  frequencies  are  controlled.  Then  a  procedure  called  interrogation  is  performed. 
The  beams,  directed  to  the  gas,  have  to  be  collinear  and  circularly  polarized.  Two 
frequency  values  are  important:  the  difference  between  the  two  beams  frequency,  and  the 
laser  center  frequency.  Both  are  selected  in  order  to  excite  the  atoms  to  an  optically 
excited  state;  when  “the  difference  frequency  between  the  two  laser  fields  equals  the 
ground  state  hyperfine  splitting,  a  coherence  is  generated  in  the  atomic  vapor”  [31]. 
Consequently,  the  atoms  that  are  coherently  excited  do  not  scatter  the  laser  beam,  and  the 
power  on  the  detector  increases  [31]. 

Two  important  technological  advances  were  necessary  to  allow  the  CSAC 
construction  and  manufacturing.  First  was  to  confine  high-pressure  cesium  vapor  at  65°C 
in  a  milimetric  cell.  Another  necessary  technological  advance  was  the  use  of  high 
bandwidth  VCSELs,  which  are  necessary  to  the  interrogation  process  [31]. 


VCSEL  A/4  Cs  Cell  Photodetector 


Figure  64.  CSAC  Using  Coherent  Population  Trapping  Technique 


Source  [31]:  R.  Lutwak,  D.  Emmons,  W.  Riley,  and  R.  M.  Garvey,  “The  Chip-Scale 
Atomic  Clock  -  Coherent  Population  Trapping  Conventional  Interrogation,”  34th  Annual 
Precise  Time  and  Time  Interval  (PTTI)  Meeting.  Beverly,  MA,  pp.  542,  2002. 
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c.  Voltage  Controlled  Oscillators 

Oscillators  are  expected  to  have  low  phase  noise,  high  Q-factor  and  small  drift; 
crystal  oscillators  are  used  because  they  meet  these  requirements.  Additionally,  some 
applications  demand  also  a  tunable  signal  reference. 

An  oscillator  that  is  designed  to  change  its  resonance  frequency  has  the  drawback 
of  being  more  susceptible  to  unintentional  tuning  [46].  Consequently,  an  oscillator  has 
stability  inversely  proportional  to  tuning  range  due  to  changes  in  resonance  element  load 
reactance;  however,  a  VCO  is  still  one  of  the  most  important  devices  for  synchronized 
systems. 

Frequency  tuning  in  a  VCO  is  performed  using  a  variable  reactance  load  on  the 
resonant  element,  which  changes  its  oscillation  period.  Usually  a  varactor  is  employed  for 
this  task.  For  relaxation  oscillators,  the  variable  capacitance  changes  the  discharge  time, 
while  LC  oscillators  change  their  natural  oscillation  frequency;  however,  there  are  many 
techniques  to  implement  a  VCO  for  different  requirements  of  phase  noise,  Q-factor  and 
accuracy. 

In  Table  14,  several  COTS  oscillators  and  their  performance  specifications  are 

listed. 
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Table  14.  Oscillator  Performance  Comparison 


Type 

Model 

Range 

(MHz) 

Phase  noise 
(dBc) 

Accuracy 

(ppm) 

Reference 

vcxo 

Crystal  CCO- 
083 

10-200 

-80  @  10  Hz 
-112  @  100  Hz 
-140  @  1  kHz 
-150  @  10  kHz 

±25 

[49] 

ocxo 

Crystal 

CO27VS05D 

E 

10 

-102  @  10  Hz 
-132  @  100  Hz 
-150  @  1  kHz 
-152  @  10  kHz 

±0.2 

[50] 

TCXO 

Crystal 

CVT32 

10-30 

-83  @  10  Hz 
-115  @  100  Hz 
-135  @  1  kHz 
-140  @  10  kHz 

±2.5 

[51] 

MCXO 

Vectron 

MX-503 

8-26 

-65  @  1  Hz 
-93  @  10  Hz 
-118@  100  Hz 
-140  @  1  kHz 
-154  @  10  kHz 
-156  @  100  kHz 

±0.5 

[52] 

VCO 

Vectron 

VS-705 

120- 

1000 

-30  @  10  Hz 
-60  @  100  Hz 
-85  @  1  kHz 
-105  @  10  kHz 
-130  @  100  kHz 

±100 

[53] 

CsFS 

Micro  semi 
5071 A 

5 

-106  @  1  Hz 
-130  @  10  Hz 
-145  @  100  Hz 
-150  @  1  kHz 
-154  @  10  kHz 
-154  @  100  kHz 

± 

0.000001 

[38] 

RbFS 

Symmetricon 

8040C 

5 

-72  @  1  Hz 
-95  @  10  Hz 
-130  @  100  Hz 
-140  @  1  kHz 
-148  @  10  kHz 

± 

0.00005 

[54] 

CSAC 

Micro  semi 
Quantum 
SA.45s  CSAC 

10 

-50  @  1  Hz 
-70  @  10  Hz 
-113  @  100  Hz 
-128  @  1  kHz 
-135  @  10  kHz 
-140  @  100  kHz 

± 

0.00005 

[55] 
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